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CHAPTER  I 


INTRODUCTION  AND  SUMMARY 


1.1  INTRODUCTION 

The  accuracy  of  advanced  ballistic  missile  systems  depends 
ultimately  on  the  capability  to  compensate  for  trajectory  errors 
encountered  in  boost,  deployment  and  reentry  phases.  The 
error  compensation  procedures,  however,  depend  on  the  extent 
to  which  system  performance  anomalies  can  be  quantified  in 
terms  of  guidance,  navigation,  or  basic  configuration  para¬ 
meters.  This  quantification  is  achieved  most  efficiently  by 
development  of  a  systematic  mathematical  modeling  methodology. 

The  reentry  phase  represents  a  mission  stage  which  is 
difficult  to  model.  This  is  because  of  the  complexity  of 
the  dynamic,  aerodynamic,  thermodynamic,  and  material 
phenomena,  all  of  which  affect  the  reentry  body  motion  and, 
hence,  accuracy.  Development  of  the  required  modeling  methodol¬ 
ogy  is  therefore  a  continuing  effort  of  some  technical  sophisti¬ 
cation.  Many  attempts,  both  theoretical  and  experimental,  have 
been  made  to  understand  the  flow  field  around  a  spinning  ballis¬ 
tic  reentry  vehicle.  Theoretical  techniques  are  difficult  to  use 
for  accurately. modeling  the  aerodynamics  because  of  the  complex 
interaction  of  the  many  physical  factors.  Wind  tunnel  tests 
and  analysis  significantly  contribute  to  understanding  these 
factors,  but  are  not  yet  able  to  provide  sufficient  data  for 
many  extreme  regimes  Ce*g*»  high  Reynolds  number,  high  windward 
meridian  rate  and  in  some  cases  high  Mach  number) .  It  is 
therefore  necessary  to  complement  such  analysis  and  wind  tunnel 
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tests  with  procedures  for  modeling  reentry  body  responses  from 
actual  flight  data.  Such  a  flight  test  data  procedure  must 
recognize  that  both  deterministic  and  stochastic  effects  are 
encountered  over  a  series  of  flight  tests.  The  procedure  must 
extract  maximum  information  from  data,  whether  it  is  from  one 
test  or  it  is  from  several  tests.  These  requirements  are  met, 
at  least  conceptually,  by  system  identification  technology. 

System  identification  is  the  process  of  extracting  mathe¬ 
matical  model  structure  and  parameters  from  a  set  of  test 
data  [1].  An  overview  of  the  essential  steps  of  this  process 
is  shown  in  Figure  1.1.  Data  is  first  conditioned  to  mini¬ 
mize  the  effect  of  noise.  Time  histories  of  unmeasured 
channels  are  then  reconstructed  from  available  measurements. 
Then,  an  algorithm  (or  series  of  algorithms)  is  used  to 
estimate  the  dynamic  and  aerodynamic  factors  which  probably 
affect  vehicle  aerodynamics.  Finally,  parameter  estimates 
are  refined.  (Details  of  this  procedure  are  given  in 
Chapter  II.) 


T 

AMcoTwt : 


Figure  1.1  System  Identification  Procedure  for 

Ballistic  Reentry  Vehicle  Data  Processing 
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This  report  summarizes  the  preliminary  application  of 
this  system  identification  methodology  to  ballistic  reentry 
vehicles.  The  ultimate  objective  of  this  application  is  the 
precise  modeling  of  the  behavior  of  this  system  to  subsequently 
systematize  its  accuracy  improvement. 

This  report  is  organized  as  follows.  Chapter  II  reviews 
relevant  equations  of  motion  of  a  ballistic  vehicle  and 
Chapter  III  discusses  the  system  identification  procedure  for 
such  vehicles  during  reentry.  Identification  results  for 
simulation  data  are  shown  in  Chapter  IV,  and  for  flight  test 
data  are  detailed  in  Chapter  V.  Chapter  VI  gives  the  conclu¬ 
sions  from,  this  effort. 

1.2  PRINCIPAL  DEVELOPMENTS  AND  RESULTS 


The  following  major  developments  and  results  have  been 
obtained  towards  the  development  of  an  integrated  system 
identification  method  to  estimate  accurate  models  for  aerody¬ 
namic  characteristics  of  ballistic  vehicles  from  limited  data: 

(1)  .4n  accurate  reconstruction  procedure  is  developed. 
This ■  combines  accelerometer  and  angular  rate  mWsure- 
ments  to  estimate  the  initial  condition  and  total 
angle-of -attack  and  windward  meridian  time  histories. 
The  in-plane  and  out-of-plane  moments  may  then  be 
estimated. 

(2)  ^  general  spline  representation  is  developed  to 
express  the  highly  transient  aerodynamics  of  the 
reentry  vehicle.  The  model  structure  determin¬ 
ation  step  is  modified  to  select  the  more  important 
terms  from  the  general  spline  representation. 

(3)  The  computer  programs  are  modified  such  that  mechan¬ 
isms  of  body-fixed  and  wind-fixed  forces  and  moments 
can  be  compared^  liT addition,  a  model  combining  the 
two  mechanisms  may  also  be  tested. 

(4)  The  nonlinear  maximum  likelihood  computer  proj^ram  is 
configured  to  work  with  general  ballistic  vehicle 
aerodynamic  models.  This  includes  body-fixed  and 
wind-fixed  forces  and  moments. 
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The  software  was  applied  to  a  simulation  data  record 
by  The  Aerospace  Corporation.  This  was  followed  by  an 
application  to  a  flight  test  data  record. 

The  following  preliminary  conclusions  are  based  on 
these  applications: 

(1)  The  algorithm  application  to  noise-free  simulated 
data  provides  basic  validation  of  the  software 
data  processing  procedure.  It  is  found  that 
relatively  simple  mathematical  models  can  effec¬ 
tively  represent  ballistic  vehicle  aerodynamic 
phenomena  for  this  simulated  trajectory  (except 
boundary  layer  transition) .  Since  the  simulated 
data  used  in  this  evalution  was  noise  free,  esti¬ 
mated  parameter  errors  are  small. 

(2)  The  algorithm  application  to  a  single  flight 
record  demonstrated  that  the  software  could  re- 
construct  ballistic  trajectory  parameters  required 
for  the  development  of  aerodynamic  models.  Aero¬ 
dynamic  models  could  then  be  estimated.  Note 
that  the  data  reconstruction  software  compensates 
for  bias  and  scale  factor  error  in  the  flight 
data  (Appendix  B) . 

(3)  The  results  of  processing  the  single  flight  record 
indicate  the  possibility  of  adequately  describing 
the  boundary  layer  transition  by  a  basically  wind- 
fixed  flow  with  a  short  time  interval  of  body- 
fixed  flow.  Definitive  conclusions  on  the  general 
nature  of  vehicle  forces  and  moments  (e.g.,  body- 
fixed  and  wind-fixed)  are  not  appropriate  on  the 
basis  of  a  single  flight  record.  Validation  of 
this  model  requires  analysis  of  further  flight 
data. 
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CHAPTER  II 


DYNAMICS  AND  AERODYNAMICS  OF  BALLISTIC  REENTRY 

•VEHICLE  FLIGHT 


2.1  INTRODUCTION 

The  complete  reentry  flight  of  a  ballistic  vehicle  is 
characterized  by  multiple  phenomena.  Above  300,000  feet,  air 
density  is  low  and  the  motion  is  almost  purely  kinematic. 

The  aerodynamic  forces  become  more  important  as  the  air  density 
increases.  The  boundary  layer  is  laminar  at  low  Reynolds  number 
flow,  then  goes  through  a  transition  regime  to  turbulent  flow 
with  increasing  Reynolds  number.  The  roll  rate  remains 
essentially  constant  during  the  high  altitude  kinematic  motion 
and  may  speed  up  or  slow  down  with  the  subsequent  aerodynamic 
forces  and  moments. 

Various  designs  of  reentry  vehicle  exhibit  significantly 
different  behavior  during  reentry.  The  behavior,  in  general, 
is  dependent  on  both  the  trajectory  parameters  (e.g.  initial 
angle-of-attack,  maximum  dynamic  pressure)  and  the  reentry 
vehicle  design  (e.g.  bluntness,  surface  roughness,  etc.).  In 
addition,  significant  trajectory  deflections  may  occur  in  the 
transition  region,  where  the  boundary  layer  changes  from  fully 
laminar  to  fully  turbulent.  To  study  these  trajectory  deflec¬ 
tions  and  related  behavior  anomalies  and  to  estimate  them  from 
dynamic  measurement  data,  a  comprehensive  set  of  vehicle  dy¬ 
namic  equations  is  required. 

The  equations  for  a  ballistic  reentry  vehicle  motion  have 
been  studied  by  several  authors  previously  [2-13].  This 
chapter  reviews  the  dynamics  of  a  ballistic  reentry  vehicle. 

The  significant  aerodynamic  forces  and  moments  are  indicated, 
so  that  the  resultant  describing  equations  can  serve  as  the 
basis  for  identifying  aerodynamic  effects  from  flight  test 
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data.  Certain  useful  approximations  to  these  equations  are 
derived  specifically  for  use  with  system  identification  tech¬ 
niques  of  later  chapters, 

2.2  DYNAMICS  EQUATIONS 

In  this  section,  we  define  reference  frames  commonly  used 
in  ballistic  reentry  vehicle  (BRV)  study  and  then  review  general 
six  degree-of-freedom  (DOF)  equations  of  motion  for  this  vehicle 
during  atmospheric  flight. 

2.2.1  Reference  Frames 


Figure  2.1  shows  three  reference  frames  which  may  be  used. 
Frame  A  is  the  wind  frame,  in  which  the  axis  is  aligned 

with  the  instantaneous  velocity  and  is  local  geocentric 

horizontal.  Frame  B  is  the  body  fixed  frame.  The  Xg  axis 


I 
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FRAME  A  -  UINO-FIXED 
FRAME  3  •  SOOT -FIXED 

FRAME  £  -  ENGINEERING 

AEROBALLISTIC  | 


I 


RELATIVE 
VELOCITT 
VECTOR,  V 


TOTAL  A.VGLI-0F.AT7ACX.  a.  IS  THE  ANGLE  BETAEEN  Xj  ANO  Xy  '..'INO'VARD 
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PLANE  ANO  Tj. 

Figure  2.1  Convenient  Frames  for  Representing  Ballistic 
Reentry  Vehicle  Aerodynamics 
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is  the  axis  of  symmetry  going  through  the  nose  and  Yg  axis  is 
aligned  with  a  specific  direction  at  some  point  in  time  and  then 
fixed  to  the  body.  The  engineering  aeroballistic  frame,  E,  is 
another  body  fixed  frame  except  that  it  does  not  roll  with  .the 
body.  The  relative  orientation  between  the  aeroballistic  frame 
and  the  body  frame  is  the  body  roll  angle.  In  the  engineering 
system,  a  yaw,  pitch,  roll  sequence  is  used  to  rotate  frame  A 
through  frame  E  to  body  fixed  frame  B.  The  three  Euler  angles 
are  0  and  (j).  In  the  classical  system  a  roll,  pitch,  roll 

sequence  is  used  to  rotate  frame  A  into  frame  B  giving  a  dif¬ 
ferent  aeroballistic  frame.  Both  frames  have  their  advantages, 
however,  we  will  only  use  the  engineering  frame  and  engineering 
Euler  angles.  The  total  angle-of-attack  a  and  the  windward 
meridian  are 

a  *  cos'^Ccos  ip  cos  0) 

.1  C2.1) 

»  <1)  -  tan  (tan  tj^/sin  0) 

The  wind  fixed  frame  is  not  inertial.  It  rotates  to  the 
ext'ent  of  trajectory  bending.  In  most  flights,  however,  the 
trajectory  bending  is  negligible  compared  to  other  motions 
and  is  often  neglected.  Unless  otherwise  mentioned,  in  all 
our  analysis,  the  wind  fixed  frame  will  be  considered  inertial. 

2.2.2  Dynamic  Equations 

The  six  DOF  equations  of  motion  of  a  ballistic  reentry 
vehicle  may  be  written  in  any  of  the  three  frames  discussed 
above.  The  wind  frame , complicates  the  kinematic  element  of 
the  dynamic  equations  but  the  aerodynamic  forces  and  moments 
are  easily  expressed  in  terms  of  components  in  the  plane  of 
the  wind  vector  and  the  vehicle  axis  of  symmetry  (called  in¬ 
plane)  and  another  component  orthogonal  to  the  plane  (called 
out-of -plane) .  This  is  useful  because  during  a  major  portion 
of  the  flight,  through  vehicle  symmetry,  the  forces  and 
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moments  are  related  to  the  direction  of  relative  velocity. 

The  on-board  measurements  are  also  difficult  to  express  in 
terms  of  relevent  dynamic  variables.  The  aeroballistic  frame 
is  also  useful  under  certain  circumstances,  for  example  it 
can  easily  depict  the  tricyclic  motion  of  a  BRV.  The  measure¬ 
ment  equations  and  the  kinematic  elements  of  the  dynamic 
equations  are  quite  straightforward  in  the  body  fixed  axis 
system.  The  aerodynamic  in-plane  and  out-of-piane  forces 
and  moments  must,  however,  be  transformed  through  time  varying 
angles  to  be  used  in  the  body  fixed  axis  system.  Chrusciel  [3] 
uses  wind  fixed  equations  of  motion.  Because  these  equations 
separate  out  the  in-plane  and  out-of-plane  aerodynamic  moments, 
they  are  easier  to  work  with  when  simple  estimation  procedures 
are  used.  However,  the  body  fixed  equations  are  more  useful 
with  advanced  identification  procedures,  because  of  the 
need  to  correlate  multivariable  measurements  as  accurately  as 
possible  and  such  accuracy  is  practically  achieved  by  using  un¬ 
transformed  measurements  in  the  body  axis  system  [1], 

The  dynamic  equations  of  motion  of  a  rigid  reentry  vehicle 
in  the  body  axis  system  are 

u  =  vr  -  wq  -  g  sin  9  +  -  C 

^  *  m  X 

V  =  wp  -.  ur  +  g  sin  <t>  cos  9  +  (2.2) 

Roo® 

w  »  uq  -  vp  +  g  cos  cos  9  +  C 

111  M 

p  ■  -r— 

X 

q  -  (1  -  I^/Dpr  +  ’(2.3) 

r  =  -  (1  -  I,^/I)pq  + 
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^  =  p  +  q  sin  (J>  tan  0  +  r  cos  (}»  tan  9 


0  *  q  cos  (J>  -  r  sin  <j)  (2.4) 

ill  »  (q  sin  4)  +  r  cos  4))/cos  0 

Note  that  <ti ,  9  and  i|i  are  body  Euler  angles  with  respect  to  an 

inertial  axis  system. 

These  equations  assume  that:  (1)  the  body  is  symmetric 

about  the  x  axis  with  moment  o£  inertia  I^,  and  (2)  the 

moments  o£  inertia  about  the  y  axis  and  the  z  axis  are 

equal  to  each  other  .  C„,  and  C  are  the  £orce  coe££icients 

X  y  z 

and  C„,  C  and  C  are  three  moment  coe££icients .  I£  the  wind- 
Z.  m  n 

£ixed  axes  are  assumed  inertial,  the  £irst  three  equations  are 
not  required  and  <j) ,  9,  lii  become  Euler  angles  with  respect  to 

the  wind  axis  system. 

u  *  V  cos  9  cos  i|; 

V  *  V  cos  ill  sin  9  cos  4)  -  V  sin  ip  sin  (p  (2.5) 

w  »  V  cos  i|i  sin  9  sin  4>  +  V  sin  ip  cos  4) 

where  V  is  total  speed.  As  mentioned  previously,  this  is  a  reason¬ 
able  assumption  in  analyzing  aerodynamic  moments  acting  on  the 
reentry  vehicle.  The  moment  coe££icients  C„  and  mav  be 

written  in  terms  o£  the  in-plane  (in  the  plane  containing  the 
wind  vector  and  the  vehicle  x-axis) ,  out-o£-plane  (plane  per¬ 
pendicular  to  the  wind  vector  and  containing  the  vehicle  x-axis), 
and  body-£ixed  moment  coe££icient.  Let  and  be 

the  wind-£ixed  in-plane  and  out-o£-plane  moments  and  body-£ixed 
moment,  respectively.  Then, 
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A 


C2.6) 


Cm  ”  C^  sin  <p^  +  C'  cos  <p^  +  cos  n 

Cn  =  C'  cos  <j>^  -  C'  sin  (J)^  +  sin  n 

where  n  is  the  relative  orientation  of  the  body-fixed  moment 
and  the  vehicle  y  axis.  The  windward  meridian  4)  ,  and  the 
total  angle-of -attack,  as  indicated  earlier,  are  functions  of 
the  three  Euler  angles.  In  addition,  9  and  ’p  are  often  small 
and  Eqs.  (2.4)  can  be  simplified  to 


•  * 

4)  »  p  +  t(; 

0  »  q  cos  ({)  -  r  sin  (f>  (2.7) 

4)  *  q  sin  <|>  +  r  cos  (p 

The  expressions  for  windward  meridian  and  total  angle-of-attack 
also  simplify  to 


(2.8) 


It  is  clear  from  the  above  discussion  that  accurate  de¬ 
scription  of  BRV  dynamic  motion  requires  models  for 
and  The  roll  equation  is  essentially  uncoupled  from  the 

pitch/yaw  equations.  Roll  rate  behavior  and  roll  rate  anomalies 
may  be  described  using  a  model  for  .  Further  approxi¬ 

mations  to  the  above  equations  may  be  made  by  assuming  a  constant 
p  (i.e.,  C£_  ■  0) .  This  reduces  the  number  of  equations  to  five 

and  makes  the  kinematic  term  in  the  q  and  r  equations  linear. 
This  is  usually  done  to  obtain  the  tricyclic  motion  of  the  re¬ 
entry  vehicle  [6]. 
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Models  for  and  are  required  to  describe  pitch/ 

yaw  motions  which  are  of  particular  interest  in  computing  the 
lateral  dispersions  of  a  BRV.  This  report  concentrates  on  tech¬ 
niques  which  may  be  used  to  develop  such  a  model  from  flight 
test  measurements.  If,  in  addition,  trajectory  bending  is 
required  to  be  determined  during  any  portion  of  the  flight, 
models  for  C„  and  C  should  also  be  identified.  The  dif- 

•  y  2 

ferences  between  the  values  of  v  and  w  obtained  by  propa¬ 
gating  Eqs.  (2.2)  and  by  solving  Eqs.  (2,5)  would  give  tra¬ 
jectory  bending  in  the  body  axis  system.  The  lateral  velocity 
components  arising  from  trajectory  bending  may  be  directly 
converted  into  miss  distance  on  impact. 

When  the  boundary  layer  is  purely  laminar  or  purely  turbu¬ 
lent,  approximate  models  for  and  are  known.  How¬ 

ever,  much  less  is  known  in  the  transition  region,  where  the 
boundary  layer  changes  from  completely  laminar  to  completely 
turbulent  over  a  period  of  time.  This  region  is  the  most 
important  from  the  viewpoint  of  the  application  of  the  state- 
of-the-art  identification  methods. 

2.3  BOUNDARY  LAYER  TR.'XNSITION 

.As  the  Reynolds  number  of  a  ballistic  reentry  vehicle  in¬ 
creases  because  of  increasing  air  density,  the  boundary  layer 
at  a  point  on  the  surface  of  the  vehicle  starts  becoming  turbu¬ 
lent.  The  transition  front  then  proceeds  to  other  parts  of 
the  vehicle  surface.  This  phenomenon  is  called  the  boundary 
layer  transition  (BLT) .  Deterministic  and  random  effects  during 
boundary  layer  transition  has  been  studied  both  theoretically 
and  experimentally  based  on  wind  tunnel  and  flight  tests  [8-15]. 

Boundary  layer  transition  front  affects  vehicle  aerodynamic 
forces  and  moments.  Transition  to  turbulent  boundary  layer  is 
accompanied  by  an  increased  skin  friction  and  an  increased 
pressure  in  the  boundary  layer.  If  the  transition  were  symmetric 


around  the  body,  no  net  moment  will  result  on  the  body,  though 
the  variation  of  in-plane  moment  with  angle-of-attack  may  change. 
In  addition,  the  out-of-plane  moment  will  Be  negligible.  The 
transition  is  rarely  symmetric,  however,  The  asymmetry  may  take 
several  forms,  one  example  of  which  is  shown  in  Figure  2.2 


Figure  2.2  Example  of  Boundary  Layer  Transition  Front 


Wind  tunnel  tests  and  other  studies  have  shown  that 
several  factors  affect  aerodynamic  flow  around  a  BRV  during  the 
BLT.  Some  of  the  more  important  variables  are: 

(a)  cone  angle  of  the  reentry  vehicle 

(b)  nose  bluntness,  ratio  of  nose  radius  to  base 
radius 

(c)  nosetip  and  heatshield  material  (this  determines 
the  extent  of  ablation  and  mass  addition  to  the 
boundary  layer) 

(d)  initial  roll  rate,  speed,  and  angle  at  reentry 

(e)  mass  unbalance 

(f)  surface  roughness  and  other  surface  nonuniform¬ 
ities  (e.g.  radar  windows,  antennas) 

Moments  on  the  vehicle  during  BLT  cause  pitch/yaw  and  roll 
rate  anomalies  which  are  not  predicted  by  inviscid  flow  theory. 
Pitching  and  yawing  motions  of  the  BRV  affect  angle-of-attack, 
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producing  nonzero  average  sideward  force  on  the  vehicle  leading 
to  trajectory  bending.  This  trajectory  bending  and  lateral 
velocity  changes  caused  by  trajectory  bending  are  two  of  the 
most  major  reasons  for  reentry- induced  impact  point  dispersions. 

The  manner  in  which  BIT  affects  vehicle  moments  depends 
upon  the  transition  pattern.  Several  models  have  been  suggested 
for  the  transition  phenomenon  [3,14],  some  of  which  are: 

(a)  unstable  damping  coefficient 

(b)  asymmetric'  in-plane  flow 

(c)  time- varying  body  fixed  moment 

(d)  total  moment  lag 

The  consequences  of  these  models  on  body  moments  are  discussed 
in  Table  2.1. 


Table  2.1 

Transition  Mechanisms  and  Their  Effects  on  Vehicle 

Aerodynamic  Moments 


tWHSITIQR  nCOWMlSH 

CrrCCT  ON  ACNOOVNAMIC  mohints 

COMCNTS 

(«)  UN!»TASiC0  OWIW 
cocFficiChr 

•  2CW  OU(-Of-n.AM  WMCNIS 

•  SCIF-EUITEO.  STDOHCir  OCPCNOCNT  ON  INITIAl  COHOITIONS 

•  NO  EFFCCr  ON  NOHCNr  vs.  TOMl.  ANCIC  OF  IfTACK  CUNVt 

•  CANNUT  LXFlAln  OUT'OF'FIANC  MOHENTS 

•  CbNTMOICTEO  IT  NANT  OISERVATIQNS 

•  LOOKS  ACASONAILE  ON  A  FIRST  GLANCE  AT  THE 
TIME  HISTORIES 

(fel  ASytMTRIC 

IN-PIAM(  flOM 

f  IN-KANf  OCPfNOCNT  ON  TOTAi  AMLC-OF-AHACK 

ANO  NCTNOIOS  MMCR 
«  2CM  (Uir-OF-PIAIIC  NQMCNrS 
•  NON’MOT  FIICO 

•  CANNOT  EIFIAIN  OUT-OF-PIAHE  MONCNTS 

•  WNCNTS  DO  NOT  DEFEND  UFON  BOOT 

ORIENTATION 

(c)  6001  IMCOIURITICS 
CAUSC  ASYWETRIC 

Finu  MTTCM 

•  UOY-FlttO  IIMC  VARVIM  WVCNT.  THC  AZIWYN  OF 

loor^Fiieo  NWCNr  mi  otMice  with  timc 

•  OUT-nF.FlANC  rUKNT  CONFONCNTS 

•  NINON  CHAMCS  IN  MMCNT  V6.  AMLC-OF-ATTACK  CUNVf 

•  CAUSCO  IT  SURFACE  ROUGHNESS.  SURFACE 
ASVWETRY,  ANTENNAS.  JOINTS.  CAUSING 
ACCELERATEO  TRANSITION  TO  lURIULENT 

BOUNOART  LATER.  DIFFERENT  IRRECUURITIES 
mi  IE  PRONINANT  AT  DIFFERENT  TINES. 

(ill  tcirw.  HKIir  IM 

•  TOTAL  ANTMlf  OCFCNOCNT  ON  ACNOOYRAMIC  VAilAMC 

•  MiC  LAO  OCFCNOCNT  ON  UifKMANO  MCRIOIAH  MTC 

•  N(M-U«  OUT-OF-FIANC  NNNCNT 

•  MTHENT  LA«  NAT  DEFEND  ON  UINOMARO  TCRIOIAN 
RATE  ANO  OTHER  VARIAILES 

Time-varying  body  fixed  moments  are  caused  by  small  or 
large  body  irregularities.  Small  body  irregularities,  like  sur¬ 
face  roughness  and  joints,  etc.,  cause  random  asymmetries  in 
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BLT  front.  Macroscopic  irregularities  like  antenna  windows  and 
known  vehicle  asymmetries  produce  deterministic  effects  on  BLT. 
Though  this  moment  is  fixed  in  the  body,  its  magnitude  may  change 
with  time.  In  general,  such  a  body  fixed  moment  will  lead  to 
both  in-plane  and  out-of -plvine  aerodynamic  components.  With  le- 
ference  to  Figure  2.3  these  components  are: 

“  ^b^^^  cos[^j,-*„(t)] 


b  ^w’ 


*  C^Ct)  sin[(<i^-<^^(t)  ] 


POINT  OF 

APPLICATION  OF  THE 
BODY -FIXED  MOMENT 


AND  ARE  THE  PROJECTIONS  OF  THE  WINO-FIXED  AXES  ON  THE 
Yg-Zg  PLANE. 


Figure  2.3  Body-Fixed  Moment  Mechanism  to 
Explain  the  BLT 


•Note  that  is  a  constant  and  {}),^(t)  is  a  function  of  time.  The 
values  of  C^(t)  would  depend  upon  the  aerodynamic  variables  (e.g. 


total  angle -of -attack  and  Reynolds  number)  and  the  relative 
orientation  of  the  location  of  the  body  fixed  moment  with  wind 
direction.  The  in-plane  moment  will  also  have  another  component 
resulting  from  wind  fixed  flow. 

The  total  moment  lag  mechanism  is  illustrated  in  Figure  2.4. 
The  lag  angle  is  a  constant  function  of  time-varying  windward 
meridian  rate  and  the  total  moment  depends  upon  the  aerodynamic 
variables : 

’“N  ■  ‘■m 

H  *  “"*1 

Cj^Ct)  -  /Cj^^(t)  +  c^^Ct) 

-1 

c;:  “ 

M 

Cj„Ct)  »  qSc  {C^[aCt-T),  RgCt-r)]} 


Figure  2.4  Total  Moment  Lag  Mechanism  for  BLT 
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2 . 4  SUMMARY 


This  section  describes,  in  brief,  the  dynamics  and  aero¬ 
dynamics  of  flow  around  a  BRV  during  reentry.  Special  attention 
is  given  to  the  BIT  regime  where  the  boundary  layer  changes 
from  fully  laminar  to  fully  turbulent.  The  BLT  regime  is  a 
combination  of  deterministic  and  stochastic  effects.  Though 
temporal  and  spatial  random  variations  in  flow  account  for  some 
of  the  behavior  observed  during  BLT,  they  are  not  useful  because 
of  their  lack  of  predictive  capability.  The  next  section  will 
describe  state-of-the-art  system  identification  methods  which 
may  be  used  to  model  deterministic  components  of  "anomalous" 
aerodynamic  moments,  which  are  not  predicted  by  inviscid  theory. 
Future  sections  will  then  discuss  the  results  of  the  application 
of  the  identification  techniques  to  simulated  and  flight  data. 
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CHAPTER  III 


SYSTEM  IDENTIFICATION  TECHNOLOGY  FOR 
BALLISTIC  REENTRY  VEHICLES 


3.1  INTRODUCTION 


This  section  provides  an  overview  of  advanced  system  iden¬ 
tification  techniques  for  estimating  aerodynamic  coefficients 
of  ballistic  reentry  vehicles  (BRV)  from  dynamic  flight  measure 
ments.  The  following  major  factors,  characteristic  to  BRV, 
dictate  requirements  and  limitations  on  applicable  system 
identification  techniques  for  the  BLT  regime: 

(1)  Moment  and  force  coefficients  are  unknown  nonlinear 
functions  of  aerodynamic  variables  (e.g.,  angle- 

of  attack  and  Reynolds' number) . 

(2)  The  form  of  dependence  between  the  force  coefficients 
and  aerodynamic  variables  is  not  known.  Several 
qualitative  hypotheses  (some  of  which  are  conflict¬ 
ing)  have  been  proposed  to  describe  the  flow  field 
under  various  conditions. 

(3)  There  is  no  control  input  to  the  system.  All  state 
excursions  (which  are  necessary  for  extraction  of 
aerodynamic  coefficients)  result  from  non-zero 
initial  conditions  or  aerodynamic  model  changes. 

(4)  Aerodynamic  phenomena  are  possibly  a  combination 
of  deterministic  and  random  effects. 

(5)  Certain  key  variables,  like  total  angle-of -attack 
and  windward  meridian,  cannot  be  measured  directly 
and  must  be  derived  from  other  measurements.  This 
causes  correlation  between  errors  in  derived  angle 
of  angle  and  other  measurements. 

(6)  The  instrument  outputs  are  contaminated  with  random 
noise,  bias  and  other  errors. 
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General  techniques  for  the  development  o£  mathematical 
models  from  test  data  have  been  discussed  previously  [1,16, 

17] .  In  the  following,  we  show  how  these  methods  are  modified 
to  be  more  responsive  to  the  requirements  given  above.  The 
modified  procedures  are  applicable  equally  to  anomalous  regions 
like  boundary  layer  transition  and  to  other  regions  where  the 
aerodynamics  are  not  as  complex. 

Figure  3.1  gives  the  functional  flow  chart  for  identifying 
BRV  aerodynamic  models.  A  significant  issue  in  application  of 
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Figure  3.1  System  Identification  Approach  for  Ballistic 
Reentry  Vehicle  .Aerodynamic  Model  Estimation 
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system  identification  to  flight  data  is  the  degree  to  which 
the  following  steps  are  simultaneously  performed: 

(1)  reconstruction  of  system  states  includes  digita.l 
filtering  of  the  data,  estimation  of  unmeasured 
states,  and  estimation  of  aerodynamic  force  and 
moment  time  histories  by  open-loop  integration; 

(2)  Model  Structure  Determination  CMSD)  by  regression 
and  hypothesis  testing  methods;  and 

■  (3)  parameter  identification  of  the  aerodynamic  models 
using  advanced  algorithms  which,  in  addition  to 
parameters  of  aerodynamic  coefficient  models,  compute 
estimates  of  aerodynamic  derivatives,  their  un¬ 
certainties,  and  sensor  errors  and  their  uncer¬ 
tainties  . 

There  is  a  procedural  ambiguity  in  performing  these  steps 
because  each  uses  theoretical  and  computational  techniques  com¬ 
mon  to  the  other  two.  It  is  possible,  for  example,  to  perform 
reconstruction  of  system  states,  model  structure  determination, 
and  parameter  identification  in  one  algorithm.  For  example,  a 
maximum  likelihood  algorithm  can  be  structured  to  perform  all 
steps  simultaneously.  Such  a  structure,  however,  may  involve 
significant  computation  burden  which  could  be  avoided  by  recog¬ 
nition  of  specific  system  characteristics  such  as  the  effect 
of  sensor  errors  relative  to  aerodynamic  uncertainties. 

For  the  ballistic  reentry  vehicle  system,  the  optimum  al¬ 
location  is  based  primarily  on  the  relatively’ high  quality  of 
on-board  sensors,  and  the  resultant  small  errors  obtainable  in 
accelerometer  and  rate  measurements  (after  the  proper  corrections 
have  been  applied) .  For  small  errors,  the  state  reconstruction 
is  optimally  performed  separately  from  the  model  and  parameter 
estimation  steps.  Therefore,  digital  filtering  techniques  com¬ 
bined  with  efficient  integration  routines  are  suitable  for  this 
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BRV  state  reconstruction  procedure.  Model  structure  determi¬ 
nation  and  parameter  identification  are  best  performed  with  a 
maximum  likelihood  identification  approach. 

Subsequent  sections  of  this  chapter  describe  briefly  the 
above  steps  specifically  for  application  to  a  ballistic  reentry- 
vehicle  aerodynamic  characteristics  estimation.  Details  of 
the  theory  on  which  these  steps  are  based  may  be  found  in 
Refs.  1  and  16. 

3.2  DATA  RECONSTRUCTION 

Data  reconstruction  is  the  process  of  computing,  from  the 
measurements,  the  unmeasured  quantities  necessary  for  the  iden¬ 
tification  stage.  In  general,  it  produces  approximate  estimates 
of  instrument  errors  and  a  smoothed  estimate  of  vehicle  velocity 
components.  Either  concurrently  or  in  a  second  stage  the  aero¬ 
dynamic  coefficients,  total  angle-of-attack  and  other  variables 
of  interest  are  evaluated  with  appropriate  kinematic  relation¬ 
ships.  The  objective  of  methods  reported  here  is  not  to  re¬ 
assess  current  techniques  of  obtaining  the  best  estimated  tra¬ 
jectory  (BET),  nor  to  integrate  the  estimation  of  the  aero¬ 
dynamics  with  the  trajectory  reconstruction  and  instrument 
error  estimation.  Rather,  the  objective  is  to  model  the  non¬ 
linear  aerodynamics  as  definitively  as  possible,  such  that  the 
model  based  on  the  reconstruction  will  help  interpret  the 
physical  phenomenon. 

Typical  dynamic  measurements  available  from  a  BRV  flight 
are:  fa)  three  axis  rate  gyros,  (b)  three  axis  linear  acceler¬ 
ometers,  and  (c)  total  velocity,  altitude,  dynamic  pressure  and 
Reynolds  number  from  off-board  measurements  and  BET.  The  recon 
struction  procedure  is  based  on  these  measurements,  but  can  some 
times  be  modified,  if  some  of  these  measurements  are  not  avail¬ 
able. 


The  reconstruction  problem  consists  of  using  available 
measurements  to  compute  the  following  time  histories:  (a)  total 
angle-of-attack,  (b)  windward  meridian  and  rate,  and  (c)  force 
and  moment  coefficients  time  histories.  This  problem  may  be  sub¬ 
divided  into  two  steps:  (1)  estimation  of  initial  conditions, 
and  (2)  estimation  of  time  histories  given  the  initial  condition. 
These  steps  may  be  attempted  sequentially  or  concurrently. 

The  reliance  of  the  reconstruction  on  the  BET  affects  the 
model  structure  determination  phase  but  not  the  parameter  esti¬ 
mation  phase  because  the  maximum  likelihood  algorithm  is  used 
with  on-board  measurements  such  that  the  instrumentation  errors 
can  be  identified  as  well  as  the' aerodynamic  coefficient  par¬ 
ameters  . 

3.2.1  Estimation  of  Intial  Conditions 

Two  methods  were  investigated  to  find  the  initial  wind  plane 
orientation.  The  first  method  estimates  the  wind  plane  orienta¬ 
tion  prior  to  boundary  layer  transition.  In  the  second  method, 
the  wind  plane  orientation  and  time  histories  of  important  vari¬ 
ables  are  jointly  estimated.  One  method  to  determine  wind  plane 
orientation  and  total  angle-of-attack  prior  to  BLT  is  given  in 
Appendix  A  and  another  method  has  been  described  by  Belknap  [22]. 


3.2.2  Estimation  of  Time  Histories 


There  are  two  sets  of  kinematic  equations  which  are  useful 
in  time  history  reconstruction.  Using  Eqs.  (2.2  and  (2.4) 
and  neglecting  gravity 
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In  the  reconstruction  procedure,  noisy  values  of  u)  .and  A 

eg 

are  substituted  in  Eqs .  (3.1)  and  (3.2).  One  of  three  methods 
described  below  may  be  used  to  obtain  the  total  angle -of -attack 
and  the  windward  meridian  time  histories  using  these  equations. 


Method  1 

The  body  fixed  system  is  aligned  with  the  body  at  the 
initial  point  such  that  the  x-axis  lies  along  the  axis  of 
symmetry  and  the  y-axis  points  along  the  input  direction  of 
y  accelerometer.  The  three  velocity  components  are  determined 
at  the  initial  point  from  the  algorithm  of  Section  3.2.1. 


Equation  C3.1)  is  then  propagated  to  obtain  the  time  histories 
of  u,  V  and  w.  Alternatively,  only  the  v  and  w  equations  may 
be  propagated  and  the  total  velocity  used  from  off-board  measure¬ 
ments.  The  total  angle-of-attack  and  windward  meridian  are 
computed  from  v  and  w  as  follows 


=■  arctan  (w,v)  (3.4) 

Windward  meridian  rate  is  obtained  by  differentiating  Equation 
(3.4)  and  substituting  from  Equation  (3.1)  where  necessary. 

If  the  inertial  axis  is  fixed  to  the  initial  direction 
of  the  velocity  vector,  then  propagation  of  Equation  (3.2) 
gives  the  orientation  of  the  BRV  in  the  inertial  axis  system. 

A  major  problem  with  this  technique  is  that  Equation  (3.1) 
are  neutrally  stable,  therefore  measurement  errors  tend  to  build 
up  in  integration.  (Note  that  the  transition  matrix  Q  is 
singular.)  The  other  two  methods  attempt  to  avoid  this 
problem  by  assuming  that  the  trajectory  bending  is  insignificant. 

Method  2 

The  trajectory  bending  during  any  short  portion  of  the 
flight  is  assumed  insignificant  and  may  be  neglected  for  the  pur¬ 
pose  of  estimating  the  aerodynamic  moment  coefficients.  Eq. 

(3.2)  may  then  be  used  with  Eqs.  (2.5)  and  Eqs.  (3.3)  and  (3.4) 
for  estimation  of  total  angle-of-attack  and  windward  meridian. 


This  is  the  angular  rate  method  of  Belknap  [22].  If  a  small 
angle  assumption  is  also  made,  Eqs.  (3.2)  are  used  with  Eq. 
(2.8)  repeated  here  for  convenience; 


a  =  /e^  +  ilj^ 

(p^  =  arctan  (ti»,0) 


(3.5) 


Method  3 

This  method  also  assumes  no  trajectory  bending  and  uses 
both  Equations  (3.1)  and  (3.2).  Equation  (2.5)  are  applied 
to  set  up  a  filtering  formulation  in  which  results  of  (3.1) 
and  (3.2)  are  cross  checked  with  each  other  and  corrected  if 
necessary  at  each  measurement  point.  This  is  the  most  accurate 
technique  and  can  effectively  filter  out  measurement  errors 
in  the  accelerometers  and  rate  gyros. 

3.2.3  Joint  Estimation  of  Initial  Conditions  and  Time  Histories 


Advanced  filtering  techniques  enable  simultaneous  estimation 
of  initial  conditions  and  time  histories  of  important  aerodynamic 
variables.  This  is  accomplished  by  using  Method  3  above  and  con¬ 
sidering  initial  conditions  as  unknowns  in  the  filter.  Because 
the  equations  are  nonlinear  and  coupled,  several  filter  passes 
through  the  data  are  necessary. 

3.2.4  Estimation  of  Force  and  Moment  Coefficients  Time  Histories 


Force  coefficients  are  evaluated  in  the  body  frame  using 
body  accelerations  at  the  center-of-gravity .  They  are  then 
rotated  to  an  in-plane  and  out-of-plane  component  with  the 


windward  meridian.  The  body  frame  moment  coefficients  are 
computed  by: 


The  moment  coefficients  C„  and  are  also  rotated  using  the 

m  n 

windward  meridian  angle,  to  an  in-plane  moment  and  out-of- 
plane  moment  C^- 

3.3  MODEL  STRUCTURE  DEVELOPMENT 

The  model  structure  development  is  an  extremely  important 
stage  in  the  overall  system  identification  approach  for  BRV 
flight  data  processing.  A  good  model,  at  minimum, 

(1)  gives  insight  into  the  applicable  physical  phenomenon,  (2) 
explains  the  anomalous  moment  (moments  which  cannot  be  explained 
by  inviscid  flow),  and  (3)  is  capable  of  predicting  the 
nonlinear  aerodynamics  in  the  neighborhood  of  the  flight  tra¬ 
jectory.  The  model  obtained  from  a  single  flight  test  can 
only  describe  those  parts  of  the  flight  regime  through  which 
the  vehicle  travelled.  Therefore,  a  model  for  BRV  flight 
phenomena  will  not  be  complete  without  extensive,  inter¬ 
active  evaluation  of  many  flights.  However,  as  the  results 
reported  in  the  following  chapters  indicate,  it  is  possible  to 
develop,  based  on  one  flight,  models  which  interpret  the  moment 
behavior.  These  results  can  also  offer  inputs  for  the  design 
and  deployment  of  future  experimental  flights, 


51 


The  model  structure  determination  problem  is  to  explain 
the  behavior  of  aerodynamic  coefficients  time  histories  in 
terras  of  independent  aerodynamic  variables  of  interest,  functions 
of  these  variables,  or  variables  which  could  explain  a  par¬ 
ticular  physical  phenomenon  (e.g.  a  nonsymmetric  flow  sticking 
to  the  BRV  surface  during  BLT)  • 

In  general,  model  structure  determination  may  be  divided 
into  three  distinct  subproblems  (see  Ref.  29  for  details): 

(1)  Selections  of  general  functional  forms  to  relate 
the  moments  and  forces  to  independent  variables 
for  use  with  noisy  test  data. 

(2)  Criteria  to  compare  competing  models  in  order  to 
determine  the  model  with  the  best  predictive  and 
other  capabilities  for  a  class  of  responses. 

(3)  Efficient  computational  techniques  for  determining 
the  models  from  test  data. 


3.3.1  Selection  of  General  Functions 


The  first  subproblem  deals  with  determining  the  least 
complex,  but  most  effective,  model  for  predicting  system 
response.  The  linear  model,  for  example,  is  often  the  least 
complex  of  models  (if  the  number  of  required  linear  terms  is 
not  excessive),  but  will  be  limited  in  predicting  nonlinear 
response  characteristics.  The  nonlinear  model,  alternately, 
may  be  more  complex,  but  has  better  predictive  capability. 
The  most  general  formulation,  therefore,  involves  selection 
of  general  nonlinear  functions  which  can  approximate  a  wide 
class  of  unknown  linear  or  nonlinear  relationships.  Taylor 
series  polynomials,  orthogonal  polynomials,  and  Chebyshev 
polynomials  have  been  used  previously  with  good  success  in 
limited  examples.  Recent  results,  presented  in  this  report, 
show  that  a  new  formulation,  based  on  splines  with  optimal 
knots,  provides  nonlinear  mo  is  of  high  generality. 
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The  advantage  o£  spline  functions  is  that  a  certain  set  . 
of  discontinuities  can  be  included  allowing  the  model  to  rep¬ 
resent,  and  hence  interpret,  rapid  changes  in  the  underlying 
physical  phenomenon.  A  representation  for  the  spline  function 
Sjjj,  V  (x)  used  to  represent  the  true  aerodynamic  coefficient 

C  (x)  is  due  to  Greville  [23] 
aero 
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where  is  the  truncated  power  function 
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(3.7) 


where  m  is  the  polynomial  order  and  the  order  of  con¬ 
tinuity. 

The  specification  of  the  approximating  splines  requires 
the  selection  of  four  variables:  (a)  m  and  v,  (b)  number  and 
position  of  knots,  (c)  specific  terms  which  could  be  deleted, 
and  (d)  coefficients  .  The  spline  representation  is  linear 
in  unknown  parameters  .  This  facilitates  applications  of 
the  techniques  discussed  in  the  following  subsections  (see 
also  Ref.  24) . 


3 .3  . 2  Model  Selection  Criteria 


Development  of  useful  criterion  to  compare  competing 
models  is  probably  the  most  difficult  of  the  three  problems. 
Several  criterion  for  application  to  BRV  aerodynamic  model 
development  were  reviewed  and  extended  in  this  work.  These 
include : 
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(a)  residual  sum  of  squares, 

(b)  adjusted  residual  sum  of  squares, 

(c)  prediction  error, 

Cd)  Mallow's  Cp, 

(e)  total  equation  F-ratio, 

Cf)  average  expected  variance,  and 
(g)  Allen's  prediction  error  at  a  point. 

All  these  criteria  have  advantages  and  disadvantages.  For 
detailed  discussions  see  Ref.  25.  Total  F-ratio  was  used 
as  the  primary  criterion  in  this  application  and  the  pre¬ 
diction  error  was  used  as  a  cross  check.  The  importance  of 
engineering  judgement  and  its  application  in  choosing  good 
subsets  as  models  is  discussed  and  illustrated  with  simu¬ 
lation  data  in  Section  IV  and  with  flight  data  in  Section  V. 


3.3.3  Computation  Techniques 

Efficient  computational  techniques  are  an  integral  part 
of  a  successful  MSD  procedure.  Often  there  are  many  data' 
points  and  many  possible  independent  variables  or  functions 
of  independent  variables  which  may  explain  the  dependent 
variable.  To  implement  the  MSD  procedure,  therefore,  the 
computation  method  must  require  reasonable  computation  time 
and  storage.  In  engineering  practice,  a  suboptimal  model  is 
often  acceptable,  if  the  computational  time  required  to  obtain 
a  slightly  better  model  is  too  high.  Two  methods  have 
been  successfully  applied,  both  based  on  the  regression  technique. 
The  first  method  takes  less  computational  time  and  storage, 
but  does  not  necessarily  give  the  optimal  model.  Stepwise,  sub¬ 
set  regression  [26]. is  an  iterative  procedure  which  starts  by 
including  the  most  highly  correlated  variable  in  the  regression. 

At  every  step,  hypothesis  testing  is  done  to  determine  if  any 
new  variable  should  be  included  in  the  regression  or  if  any  of 
the  variables  already  in  the  equation  should  be  dropped.  The 
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partial  F-value,  defined  as  the  ratio  of  improvement  in  fit 
produced  by  entering  a  variable  in  the  regression  to  the  rotal 
fit  error,  corrected  for  the  degrees  of  freedom,  is  computed 
for  each  parameter.  The  cut-off  F-value  is  determined  from 
the  F-distribution  based  on  a  certain  confidence  level.  A 
comparison  of  the  partial  F-value  with  the  cut-off  value  gives 
the  parameters  which  may  be  included  in  the  equation  or  dropped 
from  the  equation.  The  values  of  the  coefficients  are  computed 

simultaneously . 

The  optimal  regression  technique,  called  "regression  by 
leaps  and  bounds"  [27]  finds  the  best  regressions  without 
examining  all  possible  combinations  and  gives  the  best  sub¬ 
sets  for  each  number  of  parameters  included  in  the  model. 

The  subsets  are  evaluated  based  on  the  criteria  given  in 
Section  3.3.2  to  determine  the  optimal  size. 


3.4  MAXIMUM  LIKELIHOOD  TECHNIQUE  FOR  PARAMETER  IDENTIFICATION 

The  model  structure  development  phase  gives  good  models 
for  aerodynamic  coefficients  but  the  values  of  the  parameters 
are  not  accurate.  In  addition,  the  MSD  phase  uses  reconstructed 
data  which  may  be  contaminated  with  instrumentation  errors  and 
BET  estimation  errors.  The  maximum  likelihood  parameter  iden¬ 
tification  step  estimates  the  parameters  of  the  model  from  the 
previous  step.  Significant  instrument  errors  are  estimated 
simultaneously  to  minimlte  their  effect  on  parameter  estimation 

errors . 

The  maximum  likelihood  method  is  very  general  and  compre¬ 
hensive.  Conceptually,  it  car.  be  summarized  as  follows: 

"Find  the  probability  density  functions  of  the  observations 
for  all  poLible  combinations  of  unknown  parameter  values. 
Select  thfdensity  function  whose  value  is  highest  among 
an  dLsity  functions  at  the  measured  values  of  the  ob¬ 
servations.  The  corresponding  parameter  values  are  th 
maximum  likelihood  estimates. 
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Figure  3.  3  Maximum  Likelihood  Estimates 

Suppose  9  can  take  three  possible  values:  9^,  9^,  and 

9^.  Let  the  probability  density  functions  of  observations  z 
for  these  three  values  of  9  be  as  shown  in  Figure  3.3.  Then, 
if  the  actual  observation  is  z,  82  maximum  likelihood 

estimate  of  9.  In  practice,  starting  from  a  priori  estimates, 
the  parameters  are  updated  so  that  the  value  of  the  resulting 
density  function  at  the  observations  increases  monotonically 

Suppose  we  make  a  series  of  observations  y(l),  y(2) . .y(N) , 
called  Y(N),  which  are  functions  of  an  unknown  parameter  vector 
9.  Then  the  likelihood  function  for  parameter  9  is  the  con¬ 
ditional  probability  density  function  of  Y(N)  given  9.  The 
negative  log- likelihood  function  (NLLF)  ,  -log  p(YCN)/9),  may 
be  written  as 

NLLF  *  -  log  p(Y(N)[9)  ^ 

N 

»  -  log  n  p[y(i) |Y(i-l)  ,9]  (3.9) 

N 

=  -  Z  log  p[yCi) lY(i-l)  ,9] 
i  =  l 
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This  formulation  is  particularly  useful  for  dynamic  systems 
where  the  observations  yCD.  y(2)..y(n)  can  be  ordered  with 
increasing  time.  The  conditional  probabilities  lead  directly 
to  a  two-stage  formulation,  the  first  stage  of  which  is  a 
Kalman  filter,  dependent  on  unknown  parameter  values,  and’ 
the  second  stage  is  a  function  of  the  Kalman  filter  output. 

It  is  straight-forward  to  show  that  [16],  for  white  gaussian 
process  and  measurement  noise  sources 


NLLF  »  S  {v'^(i]B‘^(i)v(i)  +  log|BCi)|} 
i^l 


(3.10) 


where  \j(i)  are  the  innovations  at  measurement  point  y(i)  and 
B(i)  is  the  covariance  of  v(i).  Both  v(i}  and  B(i)  are  outputs 
of  the  Kalman  filter  and  are,  in  general,  functions  of  unknown 
parameters.  The  procedure,  based  on  maximizing  the  likelihood 
function  or  minimizing  the  NLLF,  therefore,  involves  minimizing 
a  function  of  the  output  of  a  Kalman  filter,  which  depends  upon 
unknown  parameter  values.  If  the  gradient  based  techniques  are 
used,  we  must  compute  the  first  and  possibly  the  second  gradient 
of  NLLF.  Exact  expressions  for  this  are  complex  but  several 
approximations  enable  an  engineering  solution.  First  B(i)  is 
assumed  time  invariant  and  is  estimated  by  differentiating  NLLF 
with  respect  to  B. 
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This  reduces  NLLF  to, 
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The  fir:?t  and  second  gradients  of  the  simplified  NLLF  are, 
approximately , 
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Thus  we  need  the  gradients  of  the  innovations  as  time  functions 
which  requires  the  propagation  of  the  Kalman  filter  and  its 
sensitivities.  The  gradient  method  iterations  are  continued 
until  NLLF  cannot  be  decreased  by  changing  the  parameters.- 
Ref.  30  discusses  the  computational  problem  in  detail. 

For  the  BRF  problem,  the  actual  process  noise  (gusts, 
random  atmosphere  density  changes,  etc.)  is  quite  low  such 
that  ML  estimation  with  an  output  error  formulation  should  work 
well.  The  software  adapted  for  the  BRV  problem  is  versatile. 
Aerodynamic  coefficients  ,  Cy^ ,  Czq,  ,  C',  are 

linear  or  nonlinear  functions  represented  by  a  polynomial 
expansion  of  the  form  ‘  . 
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Z  P  n  (z.  -  z.  )  ^^ 
^  i=*l  ^0 


where 

Pj  is  a  parameter  which  may  be  identified  by  the 
program 

z^  is  an  expansion  variable  (a  function  of  states), 
which  varies  from  region  to  region 

z^  is  the  reference  value  of  the  ith  expansion 
°  variable,  and 


n..  are  integer  exponents  of  the  expansion  variables. 


This  type  of  representation  of  completely  different  nonlinear 
models  specified  comp^'etely  by  input  has  been  adapted  to  include 
spline  function  models.  The  purpose  is  still  the  same,  to  refine 
the  spline  coefficients,  and  the  program  can  also  modify  the 
model  structure  somewhat  by  identifying  certain  knot  locations 
more  accurately. 

The  coefficient  time  histories  and  their  variances  are 
computed  from  the  optimal  parameter  value  of 

T  * 

Cjn  -  1  Wi 

CoviC^)  -  f^Cx]M'^f(x) 

f(x)  *  [fj^Cx),  f2(x),...,f  (x)]^ 

where  f(x)  is  a  vector  of  the  functions  included  in  this 
model.  Similarly,  the  dimensional  stability  derivatives  and 
their  covariances  are  evaluated  by 

9C 

\  OsC'  -  'IsC  i 

is  the  dimensional  stability  derivative  of  with 

respect  to  x^  (x^  is  a  state  variable  or  a  function  of 
state  variables) . 

3 . 5  SUMMARY 

In  the  identification  of  nonlinear  systems,  the  model 
structure  determination  phase  is  extremely  important  and  re¬ 
quires  the  reconstruction  of  related  but  unmeasured  aero¬ 
dynamic  variables  and  coefficients.  This  phase  develops  the 
model  and  gives  initial  estimates  of  parameters  in  this  model. 

The  maximum  likelihood  parameter  estimation  stage  uses  these 
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model  structures  but  goes  back  to  the  original  measurements  to 
obtain  asymptotically  unbiased  and  minimum  variance  estimates 
of  the  parameters  (under  the  assumption  that  the  model  structure 
is  exact) . 


CHAPTER  IV 


AERODYNAMIC  MODEL  IDENTIFICATION  FROM 
SIMULATED  DATA 


4.1  INTRODUCTION 

Application  of  the  software  based  on  the  advanced  system 
identification  procedure  of  the  previous  section  to  simulation 
data  is  necessary  not  only  to  demonstrate  capabilities  of 
techniques  but  also  to  verify  the  software.  The  simulation 
data  used  in  this  study  is  obtained  from  General  Trajectory 
Simulation  (GTS),  an  advanced  computer  program  developeo 
by  Aerospace  [28].  GTS  uses  complex  aerodynamic 
force  and  moment  coefficient  models  and  can  simulate  aerody¬ 
namic  flights  of  ballistic  vehicles  under  a  variety  o  con¬ 
ditions.  Three  time  histories  of  BRV  simulated  measurements 
were  supplied  by  Aerospace  to  SCI  CVt) .  The  -simulated  data 
is  generated  for' a  ballistic  vehicle  trajectory  at  200  samples 
ter  second.  The  measurements  are  noise  free.  The  data  tape 
‘also  contains  aerodynamic  variables  not  available  in  flight. 
These  variables  are  used  to  verify  the  software  but  are  not 
used  in  the  identification  as  such. 

The  following  sections  describe  the  results  of  processing 
a  data  record.  Time  histories  of  a  typical  angular  “ 

a  typical  linear  accelerometer  measurement  are  shown  in  Figure 
4  1  Since  the  measurements  are  noise  free,  the  data  is  only 
used  to  verify  the  software,  but  it  is  not  used  to  evaluate 
coefficient  estimation  errors. 

4.2  DATA  RECONSTRUCTION 

The  rate  integration  method  is  used  to  compute  the  pitch 
and  yaw  angles.  Figure  4.2  shows  a  cross  plot  of  reconstructed 
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pitch  and  yaw  angles.  The  vehicle  nose  is  moving  in  near  per¬ 
fect  circles.  Averaging  pitch  and  yaw  angles  over  two  cycles 
determines  the  initial  wind  plane  orientation  in  the  body  axis 
system  and  the  initial  total  angle-of-attack.  Note  that  the 
assumption  of  no  trajectory  bending  is  implicitly  made  in  the 
rate  integration  method.  The  reconstructed  pitch  and  yaw  angles 
are  compared  to  the  simulation  values  in  Figures  4.3  and  4.4. 

The  fits  are  quite  good  considering  that  pitch  and  yaw  angles 
are  obtained  through  straight  integration.  Figure  4.5  compares 
the  reconstructed  pitching  moment  coefficient  and  the  corres¬ 
ponding  simulation  values.  The  reconstructed  out-of-plane 
moment  is  very  small  and  was  neglected.  The  time  histories  of 
the'  following  variables  were  generated  in  the  reconstruction 
stage : 

Cl)  pitch  and  yaw  angle, 

(2)  total  angle-of-attack, 

(3)  windward  meridian, 

(4)  windward  meridian  rate, 

(5)  in-plane  and  out-of-plane  moment  coefficients,  and 

(6)  in-plane  and  out-of-plane  force  coefficients. 

Reconstruction  was  also  performed  by  propagating  the  three 
linear  velocity  equations  in  the  body  axis  system  [Eqs.  (2.2)]. 

A  steady  divergence  in  the  estimated  velocity  component  was  noted. 
This  observation  agrees  with  the  analysis  of  Chapter  III.  It  was 
concluded,  therefore,  that  the  rate  integration  method  is  the 
preferred  technique.  This  calculation  also  verified  the  data  re¬ 
construction  software. 

4.3  MODEL  STRUCTURE  DEVELOPMENT 

The  aerodynamic  model  on  which  the  GTS  generated  simulation 
data  is  based  is  complex  and  is  not  useful  for  identification; 

'it  is  not  identifiable  from  the  dynamic  measurements  taken  on  a 
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single  flight.  The  main  reason  for  nonidentif iability  is  that 
the  GTS  model  holds  over  the  entire  range  of  flight  conditions 
while  the  dynamic  model  is  just  one  realization  through  the 
flight  envelope.  Also,  the  form  of  the  simulation  model  is  not 
known  to  SCI  (Vt) .  It  is  known,  however,  that  the  aerodynamic 
coefficients  are  functions  of  the  ingle-of -attack ,  Reynolds 
number  and  Mach  number. 

As  pointed  out  in  the  previous  chapter,  the  optimal  sub¬ 
set  regression  method  requires  the  specification  of  a  maximal 
model.  The  maximal  model  selected  here  is  a  general  polynomial 
in  total  angle-of -attack  (a),  Reynolds  number  (Re),  and  Mach 
number  (y),  shown  in  the  first  column  of  Table  4.1. 

The  models  isolated  by  the  .regression  technique  for  the 
in-plane  moment  and  force  coefficients  during  TALO  319.33  sec 
to  TALO  322.32  sec  are  indicated  in  the  second  and  third 
columns  of  Table  4.1,  respectively.  A  blank  entry  indicates 
thqt  the  corresponding  term  is  not  included  in  the  model. 

The  table  shows  that  a  model  with  seven  terms  can  explain  about 
99.99%  variation  in  the  in-plane  moment  coefficient  and  two 
terras  can  explain  almost  all  the  variation  in  the  in-plane 
force  coefficient.  It  indicates  that  even  though  the  normal 
force  coefficient  depends  primarily  on  angle-of-attack,  the 
center  of  pressure  is  a  more  complex  function  of  angle-of- 
attack,  Reynolds  number  and  may  have  a  weak  dependence  on 
Mach  number. 

Table  4.2  shows  the  order  in  which  the  variables  enter  the 
in-plane  moment  coefficient  equation.  Total  angle-of-attack 
is  the  most  important  variable  and  enters  the  equation  first. 
Reynolds  number  is  the  next  in  importance  and  it,  together 
with  total  angle-of-attack,  explains  about  99%  variation  in 
the  in-plane  moment  coefficient.  The  next  two  terms  increase 
the  explained  variation  to  better  than  99.95%.  The  overall 
equation  F-value  is  high  showing  good  confidence  in  the 
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Table  4.1 


Nonlinear  Terms  Considered  to  Model  the  In-Plane 

Force  and  Moment 


NONLINEAR 

TERM* 

^n' 

*=2 

Constant  .... 

-0.00101 

-0.0128 

-0.069S 

-0.859 

Re . 

0.000226 

u . 

a2  . 

a3  . 

. 

Re2  . 

-0.0000103 

Re^  . 

au3  . 

aRe . 

u2 . 

0.0125 

. 

au  . 

Reu . 

a^Re  . 

a^u . 

,u2  . 

Reu2  . 

aRe2  . 

-0.00635 

uRe2 . 

a^Re . 

a^u  . 

a*Re . 

a2Re2  . 

A2  . 

Re2u2  . 

-0.0976 

aReu . 

VARIATION 

EXPUINED; 

99.99t 

100.  OOS 

*  « 

a  ■  Oy  -  a  ; 

I  •  0.0143 

Re  »  Re^  -  He  ;  He  -  1.939  (Nonnallzed  with  Re^  -  3.71  wmion) 
u  ■  Uj  -  u  :  u  *  0.983  {Normalized  with  »  23,200  fps) 
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estimated  equation.  The  overall  F-value  for  the  in-plane 
force  coefficient  is  also  high.  These  exceptionally  high 
equation  F-values  are  obtained  because  of  noise-free  data, 

Table  4.2 

Steps  in  the  Development  of  Model  Structure  for 
In-Plane  Moment  Coefficient 


*a  =  0.0143 

Re  =  1.939  (Re^  =  3.71  Mmion) 
u  =  O.PCl  (Vq  =  23,200  fps) 

4.4  MAXIMUM  LIKELIHOOD  PARAMETER  ESTIMATION 

The  model  developed  in  the  previous  section  of  the  in¬ 
plane  force  and  moment  coefficients  is  further  refined  using 
the  maximum  likelihood  technique.  Note  that  the  out-of-plane 
moment  and  force  were  so  small  that  they  are  assumed  to  be 
zero.  The  measurements  used  in  likelihood  estimation  are  from 
the  three-axis  rate  gyros  and  linear  accelerometers. 
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Reconstructed  variables  are  not  used  in  the  maximun  likelihood 
parameter  estimation  stage  since  these  variables  are  derived 
from  rate  gyro  and  linear  accelerometer  measurements.  Pitch 
rate,  yaw  rate  and  three  Euler  angle  equations  are  integrated, 
while  roll  rate  is  taken  directly  from  the  measurement.  In 
addition,  normal  and  side  velocities  equations  are  also  inte¬ 
grated.  The  sensor  equations  currently  available  in  the  Non- 
Linear  Systems  Control  Maximum  Likelihood  Identification 
(NLSCIDNT)  program  for  use  on  ballistic  aerodynamics  are  given 
in  Appendix  B.  Any  subset  of  these  sensors  may  be  chosen  in 
any  application. 

Maximum  likelihood  parameter  estimates  are  compared  with 
regression  estimates  in  Table  4.3.  There  are  some  differences. 
The  reason  is  as  follows.  The  regression  program  used  for 


model  development  is  efficient  computationally  and  it  is  a 

> 

suboptimal  stepwise  algorithm.  This  is  confirmed  by  comparing 
the  maximum  likelihood  F- values  in  Table  4.3  with  the  order 
the  terms  entered  the  model  in  Table  4.2.  The  terms  of  highest 
F-value  (in  order  of  decreasing  F-value)  are  constant  Re,  a, 
and  oRe,  yet  the  stepwise  program  did  not  select  the  model 
in  this  order.  These  are  the  types  of  terms  thai  would  be 
expected  to  represent  the  aerodynamics  in  this  region. 

Figure  4.7  shows  the  gyro  and  accelerometer  measurements 
and  estimates  based  on  the  maximum  likelihood  model  of  Table 
4.3.  Figure  4.8  indicates  variation  in  normalized  velocity  and 
normalized  Reynolds  number  (nromalizad  to  the  first  point  on 
the  trajectory,  where  velocity  is  23,200  fps  and  Reynolds  number 
is  3.71  million) . 

Figure  4.9  shows  the  true  value  of  the  in-plane  moment 
and  Figures  4.10  through  4.12  show  the  estimated  in-plane 
moment  and  its  stability  derivatives  with  respect  to  the  total 
angle-of-attack  and  Reynolds  number.  Also  shown  are  two 
standard  deviations  of  estimation  errors.  The  uncertainty 
in  the  stability  derivatives  represents  the  uncertainty  in 
estimating  away  from  the  simulation  trajectory.  This  uncer¬ 
tainty  quantifies  the  predictive  character  of  the  model.  Both 
stability  derivatives  of  Figures  4.11  and  4.12  show  an  increase 
in  uncertainty  in  approximately  the  last  0.5  sec.  This  is 

deafly  the  result  of  including  the  higher  order  terms  with 

2  2  2  2 
Reynolds  number  (see  Table  4.2),  Re  ,  Re  u  ,  and  -Re  . 

In  this  latter  portion  of  the  trajectory  (t  *  321.82  -  322.32), 

2 

Reynolds  number  squared  nearly  doubles,  making  the  Re  terms 
much  more  significant.  Thus,  in  Figure  4.10,  the  total  in¬ 
plane  moment,  the  uncertainty  is  greater  during  the  latter 
portion.  If  the  data  were  extended  beyond  322.32,  the  uncer¬ 
tainties  during  the  latter  portion  of  time  would  decrease. 
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Figure  4.10  Estimate  of  In-Plane  Moment  and  2o  Bounds 


Figures  4.13  and  4.14  show  the  estimated  in-piane  force 
and  its  stability  derivative  with  respect  to  angle-o£-attack. 

The  results  are  similar  as  for  the  aerodynamic  moment  coeffi¬ 
cient. 

The  model  structure  (Tables  4.1  and  4.2)  indicates  that 
Reynolds  number  is  quite  important  in  explaining  the  variation 
of  the  in-plane  moment,  but  what  the  stability  derivatives 
(Figures  4.11  and  4.12)  indicate  is  that  the  functional  repre¬ 
sentation  of  Reynolds  number  effects  is  not  adequate.  Reynolds 
number  is  changing  gradually,  almost  linearly,  hence  Reynolds 
number  effects  can  probably  be  best  modeled  with  a  combination 
of  regular  polynomials  and  spline  functions.  When  knots  are 
taken  close  enough  together,  linear  splines  in  time  can  be 
used  to  represent  Reynolds  number  and  Mach  number  effects 
(see  Chapter  V,  Flight  Data  Processing). 

4.5  SUMMARY 

This  chapter  described  an  application  of  the  integration 
system  identification  technology  to  a  BRV  simulation  data. 

The  application  indicates  briefly  how  this  technology  may  be 
us-eful  in  isolating  aerodynamic  characteristics  of  BRV  from  a 
limited  number  of  measurements  and  also  verifies  the  software 
which  implements  the  algorithms.  Complex  simulation 
aerodynamics  can  be  modeled  adequately  by  fairly  simple  func¬ 
tions.  This  chapter  does  not  attempt  to  develop  theories  of 
the  boundary  layer  transition  phenomenon.  Such  an  objective 
requires  flight  test  data  taken  on  the  vehicles  of  interest. 
Analysis  of  the  data  from  one  particular  vehicle  for  the  specific 
purpose  of  developing  a  model  of  vehicle  behavior  during  boundary 
layer  transition  is  the  subject  of  the  following  chapter. 
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CHAPTER  V 


FLIGHT  DATA  SYSTEM  IDENTIFICATION  RESULTS 


5.1  INTRODUCTION 

Previous  chapters  described  ballistic  reentry  vehicle 
dynamic  equations  and  showed  how  system  identification  tech¬ 
niques  may  be  applied  to  estimate  various  unknown  nonlinear 
force  and  moment  coefficients  in  these  equations  based  on 
limited  time  history  data.  This  chapter  describes  the  appli¬ 
cation  of  the  techniques  to  a  particular  flight  data. 

The  vehicle  under  consideration  is  a  relatively  blunt 
vehicle  with  a  small  cone  angle.  Reentry  begins  at  a  small 
angle-of-attack  of  3“.  Boundary  layer  transition  starts 
at  90,000  ft  altitude  and  is  over  at  about  70,000  ft.  On¬ 
board  measurements  included  19  temperature,  1  acoustic,  2 
vibration,  and  3  axis  angular  rates  and  linear  accelerometer 
measurements.  In  addition,  a  magnetometer  measures  roll  rate. 
.■Mtitude  and  total  velocity  are  obtained  either  from  off-board 
radar  and  optical  measurements  or  from  Kalman  filter  estimates 
of  BET.  Time  history  flight  data  at  reentry  is  provided  by 
SAMSO/Aerospace  (ABRES)  to  SCI  (Vt) .  The  variables  on  the 
data  tape  are:  (a)  time,  (b)  altitude,  (c)  velocity,  (d)  dy¬ 
namic  pressure,  (e)  rate  gyro  (p,q,r,),  (f)  accelerometer 

(a  ,  a  ,  a  ),  and  (g)  magnetometer  roll  rate. 

A  z 

Figures  5.1  and  5.2  show  time  histories  of  the  acceler¬ 
ometer  and  rate  gyro  measurements.  The  measurements  are  noisy. 
The  boundary  layer  transition  starts  after  320  sec.  The  oscil¬ 
lations  produced  during  the  BLT  period  damp  out  in  about  4  sec. 
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5.2  DATA  RECONSTRUCTION 


The  basic  measurements  used  for  identification  are  the 
body  angular  rates  along  the  three  principal  axes  (measured 
by  gyros  and  magnetometer)  and  the  body  accelerations  along 
the  principal  axes  (measured  by  three  body-fixed  accelerometers). 
In  addition,  ground-based  radars  provide  measurements  of  total 
speed  and  altitude.  Data  reconstruction  procedures  are  des¬ 
cribed  in  Chapter  III.  It  should  be  mentioned  again  that 
the  reconstructed  time  histories  are  noisy  and  are  used  only 
in  the  model  structure  development  stage.  The  maximum  likeli¬ 
hood  method  develops  filtered  estimates  of  the  time  histories 
of  these  aerodynamic  variables. 

The  aerodynamic  variables  of  primary  interest  are: 

(a)  total  angle-of-attack; 

(b)  windward  meridian; 

(c)  windward  meridian  rate; 

(d)  total  in-plane  aerodynamic  moment; 

(e)  total  out-of -plane  aerodynamic  moment;  and 

(f)  Reynolds  number. 

The  ijj  and  9  are  reconstructed  starting  from  a  zero  value 
for  ^  and  9  at  time  318.2  and  propagating  the  kinematic 
angle  equations  (see  Figures  5.3  and  5.4).  The  cross  plot  of 
9  and  is  shown  in  Figure  5.5  (note  the  difference  in  scale 
on  the  X  and  y  axes).  Initially,  the  d/\p  cross  plots 
are  close  to  circular.  The  windward  meridian  is  estimated  at 
time  318.875  using  several  methods  (described  in  Chapter  III) 
all  of  which  give  essentially  the  same  results  (see  Chapter  III). 
It  was  therefore  concluded  that  the  windward  meridian  estimate 
is  quite  good  at  the  starting  point. 

Using  this  initial  condition,  Belknap's  [22]  method  was 
used  for  reconstruction  leading  to  the  time  histories  shown  in 
Figures  5.6  and  5.7.  The  reconstructed  variables  agree  quite 
closely  with  those  derived  by  AVCO  [14].  Reconstructed  in-plane 
and  out-of-plane  moments  are  shown  in  Figures  5.3  and  5.9. 
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Figure  5.7  Reconstructed  Total  .Angle-of-Attack 
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;ECONSTRUCTEO  TInE  HISTORIES 


S.3  MODEL  STRUCTURE  DEVELOPMENT 


Development  of  physically  meaningful  model  structures 
from  the  ballistic  vehicle  flight  data  is  difficult  because  of 
the  absence  of  good  a  priori  information.  AVCO  studies [14]  divide 
the  in-plane  moment  into  three  parts:  (1)  inviscid  moment, 

(2)  a  jump  at  the  onset  of  3LT,  and  (3)  an  unexplained  part. 

The  unexplained  part  is  quite  large  compared  to  the  total  in¬ 
plane  moment.  The  out-of-plane  moment  was  not  modeled  at  all. 

The  following  procedure  is  used  for  this  effort.  At 
boundary  layer  transition  and  just  prior  to  and  following 
boundary  layer  transition,  significant  changes  in  moments 
occur  upon  the  ballistic  reentry  vehicle  resulting  from 
continuous  or  discontinuous  changes  in  the  aerodynamic  para¬ 
meters  (e.g.  the  slope  of  the  pitching  moment  vs.  angle-of- 
attack  curve)  .  The  approach  used  here  is  to  apply  statistical 
methods  to  determine  the  effects  which  are  most  significant  in 
explaining  the  observed  ballistic  reentry  vehicle  time  histories. 
These  statistical  methods  are  implemented  in  two  model  struc¬ 
ture  estimation  programs,  the  optimal  subset  regression  and 
globally  optimum  regression. 

As  a  first  step  in  studying  the  behavior  of  the  ballistic 
vehicle,  the  in-plane  moment  from  318.75  sec  to  322.0  sec  is 
considered.  Four  possible  variations  are  allowed: 


(1) 

a 

sudden  jump 

in 

r  ' 
N’ 

(2) 

a 

ramp  change 

in 

r  . 

"N’ 

(3) 

a 

sudden  jump 

in 

(4) 

a 

ramp  change 

in 

These  variations  are  allowed  at  intervals  of  0.23  sec.  Using 
the  optimal  subset  regression  technique,  the  following  model 
is  obtained: 
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(1) 

C.,  is  constant  at  -0.527 

rad 

-1, 

9 

(2) 

there  is  a  ramp  change  in 

Cn 

of 

0,000228 

sec 

C3) 

there  is  a  sudden  jump  in 
320.75;  and 

‘-N 

of 

0.000177 

at 

(4) 

there  is  a  sudden  jump  in 

321.25. 

"n 

of 

-0.000170 

at 

The  details  of  the  model  are  given  in  Table  5.4.  This  simple 
model  explains  better  than  961  variation  in  Cj^.  Figure  ...10 

Table  5.4 

Preliminary  Model  for  Cj^  (318.875  -  322.0) 


COEF-ICIENT 

VALUE 

F-RATIO 

Constant  term 

-0.000565 

-- 

'^No 

-0.0527 

1223 

Constant  ramp  in  C^^ 

0.000228 

627 

Step  in  C|^  at  320.75 

0.000177 

188 

Step  in  C|^  at  321.25 

-0*000170 

129 

Variation  explained:  96.3^ 

Overall  F-ratio:  4000 
Mean  square  .165E-06 

Mean  square  prediction  error  =  .62E-08 
a  =  a  -  0.24" 


shows  the  reconstructed  and  predicted  by  the  model 

of  Table  5.4.  This  is  an  excellent  fit  considering  that  there 
are  significant  errors  during  the  reconstruction  stage.  The 
parameter  estimates  will  improve  in  the  maximum  likelihood 
stage  where  the  equations  are  not  integrated  open-loop. 
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The  above  calculation  indicates  that  the  in-plane  moment 
acting  at  the  vehicle  c.g.  may  be  related  to  certain  aerody¬ 
namic  variables  quite  satisfactorily.  This  provides  a  test  of 
the  algorithms  on  flight  data  and  also  shows  that  simple 
mathematical  models  may  explain  the  seemingly  complex  behavior. 
Model  structure  estimation  techniques  are  now  applied  to  study 
two  major  theories  of  BIT. 

5.3.1  Wind-Fixed  Flow  Sticking  to  the  Body  for  Short  Periods 
of  Time 

One  of  the  hypotheses  proposed  for  BIT  is  that  the  flow 
field  around  the  ballistic  reentry  vehicle  is  basically  wind- 
fixed.  However,  once  a  wind-fixed  flow  is  established,  it 
adheres  to  the  body  for  a  certain  period  of  time.  This  occurs 
because  at  positive  angle-of-attack,  the  transition  front  is 
asymmetric  between  the  windward  and  the  leeward  sides  of  the 
vehicle;  the  transition  is  more  advanced  on  the  leeward  side. 

If  the  windward  meridian  changes  fast  enough,  however,  the 
flow  cannot  re-establish  instantaneously  because  of  inertia, 
leading  to  asymmetries  about  the  wind  plane.  This  asymmetry 
would,  then,  explain  the  out-of-plane  moment. 

As  shown  in  Chapter  II,  this  mechanism  converts  the  in¬ 
plane  and  out-of-plane  moment  components  into  two  unknown 
functions  Cm  (t)  and  e(T).  This  model  would  require  relat- 
ing  Cjn^^Ct)  with  angle-of-attack  and  Reynolds  number  (or 
time)  and  finding  the  function  eCx),  which  may  be  dependent 
on  t.  This  is  a  very  complex  estimation  problem.  To  test 
the  basis  of  this  theory  of  BLT,  however,  the  above  formulation 
was  simplified.  The  in-plane  moment  was  set  equal  to 
and  the  out-of-plane  moment  was  modeled  by  a  sum  of  discrete 
delays ,  i . e.  , 
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C5.1) 


cjct)  ^  c^^(t) 

^max 


In  addition,  the  time  interval  is  divided  into  several  regions 
and  e(t,T^)  is  assumed  constant  in  each  region,  i.e., 

£Ct,T^)  »  t.  (t^) 


The  problem  now  consists  of  two  steps.  In  the  first  step,  the 
in-plane  moment  is  modeled  using  the  total  angle-of-attack  and 
Reynolds  number  as  independent  variables.  In  the  second  step, 
the  out-of-plane  moment  is  related  to  lagged  in-plane  moment 
and  the  non-zero  coefficients  (fi)  required  to  explain  the 
out-of-plane  moment  are  determined. 

Model  for  the  In-Plane  Moment 

The  in-plane  moment  is  primarily  a  function  of  angle-of- 
attack;  however,  the  function  may  change  with  Reynolds  number. 
In  addition,  sudden  changes  in  the  aerodynamic  moment  may 
result  because  of  rapid  changes  in  transition  front  over  a 
portion  of  the  body.  All  these  effects  are  included  in  the 
maximal  model  specified  for  the  regression  program.  The 
regression,  program  chooses  a  subset  of  the  maximal  model  which 
best  explains  the  reconstructed  moment  using  the  fewest  degrees 
of  freedom.  The  subset  is  further  refined  using  the  globally 
optimum  regression  program. 

Splines  are  used  to  model  possible  changes  in  (the 

portion  of  in-plane  moment  independent  of  angle-of-attack)  and 
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^  i  ^j+1 


j  *  1,2,  ...  ,j 


(5.2) 


max 


in  Cm^  (the  total  angle-o£-attack  dependent  portion  of  the 

in-plane  moment).  To  model  the  in-plane  moment  from  318.875 

sec  to  323  sec,  15  knots  are  allowed  at  318.875,  320.5,  320.8, 

320.9,  321.0,  321.1,  321.2,  321.3,  321.4,  321.5,  321.75,  322.0, 

322.25,  322.50  and  322.75.  The  maximal  model  then  allows  for 

jump  and  ramp  changes  in  Caiq  and  at  these  locations. 

The  optimal  models  for  2  through  10  variables  in  the  equation 

are  shown  in  Table  5.5.  Just  two  terms,  constant  and 

Cni|^  can  explain  69.3%  of  the  variation  in  the  in-plane  moment 

and  the  addition  of  a  jump  in  at  320.8  increases  the 

explained  variation  to  90.0%.  Going  from  4  to  5  variables, 

the  jump  in  at  320.8  is  replaced  by  a  ramp  in  Cm^  at 

318.875  and  a  jump  in  at  320.8.  The  procedure  thereafter 

adds  jumps  in  at  various  time  points  until  there  are  nine 

variables  in  the  equation.  The  tenth  variable  is  a  jump  in 

C  at  322.75. 

"'a 

Figures  5.12  and  5.13  compare  0^,^  estimates  (as  time 
functions)  for  different  numbers  of  variables  in  the  equation. 

The  estimate  changes  significantly  until  there  are  five  variables 
in  the  equation.  Beyond  that  new  terms  fine  tune  the  model.  A 
similar  behavior  is  observed  for  estimate  (Figures  5.14 

and  5.15).  It  appears  that  models  with  fewer  than  five  variables 
would  not  be  physically  acceptable.  Models  with  7,  8  or  11 
parameters  are  considered  most  useful.  A  comparison  of  recon¬ 
structed  in-plane  moment  and  in-plane  moment  based  on  a  seven- 
parameter  model  is  shown  in  Figure  5.16.  The  error  is  also 
shown.  Inclusion  of  more  terms  will  improve  the  fit  between 
the  reconstructed  time  history  and  the  model  time  history. 

Figure  5.17  shows  the  error  about  the  inviscid  moment  versus  angle 
of-attack,  and  Figure  5.18  shows  the  total  estimated  in-plane 
moment  versus  angle-of -attack. 
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Table  5.5 

Models  for  the  In-Plane  Moment 


.0002 


.0001 


C  0.0 
"o 


-.0001 


-.0002 


Figure  5.12  Estimated  C  for  2,  4,  and  5  Terms  in  Regression 

"*o 

Equation  (Significant  differences  in  these  curves 
indicates  that  models  with  2  or  4  terms  cannot 
be  believed.) 


.0002 


.0001 
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-.0001 


-.0002 


Figure  5.13  Estimated  C  for  5,  6,  and  8  Terms  in  the  Equa- 


m 


tion  (Models  with  5,  6,  or  8  terms  seem  reasonable.) 
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VARIATION  WITH  TIME 
FOR  REGRESSION  MODEL 
WITH  S  TERMS 


Figure  5.14  Estimated  for  2,  4,  and  5  Terms  in  the  Equation 

(Significant  differences  in  these  curves  indicates 
that  models  with  2  or  4  terms  cannot  be  believed.) 


•  a  VARIATION  FOR 
REGRESSION  MODEL 
^ITH  8  TERMS 


Figure  5.15  Estimated  for  5,  6,  and  8  Terms  in  the  Equation 
(Models  with  5,  6,  or  8  terms  seem  reasonable.) 


0.0005 


Table  5.6  shows  the  F-value  of  various  terms  in  models  of 
different  sizes.  The  F-value  indicates  the  relative  importance 
of  a  term  in  the  regression  equation.  ,  ramp  in  at 

318.875  ,  and  jumps  in  CniQ  at  321.8  and  322.25  are  the  most 
important  variables.  F-values  also  indicate  accuracy  with 
which  these  variables  are  estimated.  (The  standard  deviation 
of  estimating  error  on  any  parameter  is  related  to  its  F-ratio, 
as  (parameter  value)/»^arameter  F-ratio.) 


Table  5.6 

Relative  Importance  of  Terms  in  the  In*Plane  Model 


COEFFICIENT 

VALUE 

STANDARD 

DEVIATION 

F-VALUE 

-0.000972 

4110. 

0 

a 

-0.0890 

0.00139 

779. 

Ramp  in  Cj^  from  318.875 

ot 

0.0138 

0.000494 

233. 

Jump  in  C^i  at  320.8 

0.000178 

0.0000116 

521. 

Ramp  in  C^j  at  321.5 

-0.000541- 

0.0000237 

197. 

Jump  in  C|^  at  321.75 

0.000245 

0.0000174 

103. 

Ramp  in  C^  from  322.75 

0.00125 

0.000124 

97. 

Model  for  the  Out-of-Plane  Moment 

Reconstructed  out-of-plane  moment  is  shown  in  Figure  5.19. 
.A.S  explained  before,  Eq.  (5.1)  is  used  to  model  the  out-of¬ 
plane  moment.  Again,  reconstructed  moment  from  318.875  to 
323.0  sec  is  considered  and  regression  methods  are  used  to 


select  the  best  model  with  various  numbers  of  parameters  (see 
Table  5.7).  The  F-ratios  of  the  parameters  are  shown  in  the 
brackets  and  the  overall  variation  explained  is  indicated  in 
the  lower  section  of  the  table. 


Table  5.7 

Model  for  the  Out-of-Plane  Moment 


MM»  OF  FMMCTEU  IN  HOOCL 

2 

3 

4 

s 

6 

7 

a 

CoAsUnt  T«ni 

-O.OOOOOlt 

-0.0000852 

•0.0000930 

-0.0000949 

-0.0000974 

-0.0000987 

-0.0000102 

0«l<y  al  0.2  IK  it  120.75 

0. 156 
(751 

0.627 

(287) 

0.479 

(1181 

0.480 

(142) 

0.481 

(149) 

0.511 

(171) 

0.160 

(12) 

a»Ui  <H  0.2  IK  it  122.5 

>0.590 

(1951 

-1.08 

(201) 

-t.06 

(202) 

•0.879 

(124) 

•0.880 

(127) 

-0.816 

(146) 

OaliX  at  0.2  iK  it  122.25 

0.640 

(59) 

0.  S04 
(33) 

0.340 

(14) 

0.298 

(11) 

Oaliy  Of  0.125  IK  it  122.25 

0.158 

(20) 

0.149 

(56) 

0.642 

(69) 

0.048 

(148) 

Dally  or  0.125  IK  it  122.75 

-0.278 

(16) 

-0.278 

(18) 

•0.201 

(18) 

Oaliy  of  0.125  iK  it  121.75 

-0.291 

(22) 

•0.137 

(10) 

Ocliy  of  0.2  IK  it  121.0 

0.284 

(12) 

Ocliy  of  0.125  IK  it  121.5 

-0.531 

(42.2) 

VUIATION  CinAIKCO; 

8.31 

25.  n 

10.91 

12.51 

35.41 

17.11 

40. 11 

The  model  with  three  parameters  has  a  constant  term  and 
two  delay  terms.  The  model  implies  that  0.627  fraction  of  the 
flow  remains  attached  to  the  body  of  0.2  sec  starting  at  320.75 
TALO.  At  322.5,  0.590  fraction  of  the  flow  stops  remaining 
fixed  to  the  body  for  0.2  sec,  leaving  a  fraction  of 
0.627-0.590  *  0.037,  which  still  remains  fixed  to  the  body 
beyond  322.5.  The  one  standard  deviations  on  the  estimated 
0.520  and  0.590  values  (using  F-values  of  287  and  195,  respec¬ 
tively)  are  0.037  and  0.042,  respectively.  The  difference 
between  the  two  values  is,  therefore,  insignificant.  We  would 
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conclude  form  this  model  that  the  out-o£-plane  moment  can  be 
explained  by  a  constant  term  and  about  60%  wind-fixed  flow 
remaining  attached  to  the  body  for  about  0.2  sec  between  320.75 
and  322.5.  As  the  model  order  is  increased,  the  two  prominant 
delays  are  0.2  sec  and  0.125  sec.  The  overall  variation  ex¬ 
plained  with  8  parameters  is  only  40.1%.  This  is  because  of 
the  high  level  of  noise  in  the  reconstruction  in-plane  moment, 
out-of-plane  moment  and  windward  meridian. 

The  reconstructed  out -of -plane  moment  is  compared  with  the 
predicted  out-of-plane  moment  based  on  the  eight-parameter  model 
in  Figure  5.19.  The  error  is  quite  random  except  it  seems  to 
be  correlated  with  the  angle-of-attack.  This  correlation  could 
exist,  for  example,  if  the  initial  estimate  of  windward  meridian 
were  incorrect.  Including  the  angle-of-attack  in  the  regression 


Figure  5.19  Reconstructed  and  Estimated  Out- 

of-Plane  Moment 


variables  gave  results  shown  in  Table  5.8.  The  best  models 
with  two  or  three  parameters  are  the  same  as  before.  The 
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seven-parameter  model  can  now  explain  a  SI. 81  variation  in  the 
out-of-plane  moment.  The  overall  nature  of  the  model  is  the 
same,  however,  and  the  delay  of  0.125  has  been  replaced  by 
a  delay  of  0.15.  The  reconstructed  and  predicted  out-of-plfne 
moments  based  on  the  seven-parameter  model  of  Table  5.8,  are 
shown  in  Figure  5.20.  The  error  looks  random  except  for  a  few 
large  peaks,  caused  probably  by  poor  data. 


Table  5.8 

Model  for  the  Out-of-Plane  Model  with  Angle-of- 
Attack  in  Regression  Equation 


5.3.2  Body-Fixed  Irregularities  Causing  Asymmetric  Flow 

Let  Cin^(t)  be  the  aerodynamic  moment  caused  by  non- zero 
angle-of-attack.  Suppose  a  body  asymmetry  at  angle  (|)|^  from 
the  behicle  y  axis  produces  a  body-fixed  moment  C^^Ct).  Then 
the  in-plane  and  the  out-of-plane  moments  experienced  by  the 
vehicle  are 

+  C^(t)  cos[(|>^j-4)^(t)] 

Cj(t)  -  C^j(t)  sin[(j>^-(t)^(t)  ] 

The  development  of  the  model  requires  relating  C^^Ct)  to 
aerodynamic  variables  Ce*g*»  angle-of-attack  and  Reynolds 
number)  and  finding  C^(t)  and  <J>^.  It  will  also  be  useful 
to  relate  C^(t)  to  some  aerodynamic  variables.  Note  that 
(J)^  may  change  with  time,  though  it  should  be  constant  over 
small  intervals  of  time. 

This  problem  is  attempted  in  two  steps.  In  the  first 
step,  eliminated  from  the  above  equations  to  give 

CfJCt)  »  Cm^/t)  +  CjCt)  cotan  [^^^-*|)^(t)] 


Regression  techniques  are  used  to  estimate  a  model  for  Cni^(t) 
and  to  estimate  Cj^Ct)  may  then  be  outained  in  the  second 

step  as 


C^Ct)  -  /[CJ.Ct)-C„„(t)l^  . 

A 

where  ' '' '  denotes  estimates.  The  correlation  between 

and  aerodynamic  variables  could  be  checked  over  the  time  history 
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if  C^(t)  is  estimated.  Regression  methods  used  in  model  struc¬ 
ture  estimation  indicate  that  this  mechanism  is  highly  unlikely. 

Our  analysis  of  body  fixed  moment  causing  dispersions 
during  boundary  layer  transition  is  based  on  the  assumption 
that  the  point  on  the  body  where  the  body  fixed  moment  acts 
is  fixed  throughout  the  time.  This  contrasts  with  Cruschiel's 
analysis  which  assumed  that  both  magnitude  and  location  of  the 
body  fixed  force  is  constant  only  over  one  cycle.  Using  our 
assumption,  the  regression  analysis  indicates . that  a  completely 
body  fixed  moment  is  unlikely  to  explain  the  out-of-plane 
observed  moments.  This  does  not  invalidate  Cruschiel's  results 
which  used  different  assumptions.  Future  analysis  should  try 
to  correlate  the  observations  using  Cruschiel's  approximations. 

5 . 4  SU^^MARY 

This  chapter  discusses  the  application  of  advanced  system 
identification  techniques  to  a  set  of  flight  data.  The  boundary 
layer  transition  region  with  its  anomalous  moments  is  analyzed. 

It  appears  that  a  wind-fixed  flow  sticking  to  the  body  for 
approximately  0.2  sec  describes  the  BIT  observed  moments 
satisfactorily . 


CHAi^TER  VI 
SUMMARY 

6 . 1  SUMMARY 

This  report  develops  a  systematic^ procedure ,  based  on 
system  identification  technology,  for  estimating  mathematical 
models  from  ballistic  reentry  vehicle  flight  data,  Accurate 
models  resulting  from  such  analysis  will  play  a  role  in  im¬ 
proving  the  overall  accuracy  of  ballistic  missile  systems. 

The  system  identification  procedure,  developed  here,  is 
general  and,  at  least  theoretically,  extracts  all  available 
information  from  a  set  of  data.  The  reconstruction  step  pro¬ 
vides  a  reduced  sensitivity  to  information  limitation  imposed 
by  a  certain  fixed  set  of  instruments,  although  it  is  recognized 
that  additional  instruments  would  provide  more  detailed  models 
and  enhance  estimation  accuracy.  Any  a  priori  information 
could  also  be  used  to  advantage  (none  was  available  for  this 
project).  The  model  structure  development  step  provides  a 
means  to  test  a  large  class  of  models’.  Therefore,  this  step 
is  extremely  useful  in  ballistic  reentry  vehicle  flight  data 
analysis,  where  a  priori  knowledge  may  be  poor  and  one  must 
resort  to  some  amount  of  trial  and  error  in  the  development  of 
useful  models.  The  model  structure  development  phase  is  followed 
by  the  relatively  sophisticated  maximum  likelihood  parameter 
estimation  stage.  Here  the  model  is  fine  tuned  and  the  unknown 
parameters  describing  this  model  are  accurately  estimated, 
while  taking  into  account  possible  errors  in  instruments.  This 
overall  structure  divides  the  development  of  BRV  mathematical 
models  from  flight  test  data  into  a  series  of  simpler  problems. 

In  addition,  it  provides  a  flexible  mechanism  to  handle  a 
variety  of  situations  in  BRV  flight  data  processing. 

This  structure  is  used  to  process  a  flight  test  data 
record  and  a  simulation  data  record  based  on  an  advanced 
digital  simulation  [28].  Both  simulation  and  flight  test 
data  records  can  be  explained  by  relatively  simple  models. 
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Even  such  simple  models  are  often  hard  to  isolate  using  clas¬ 
sical  means  because  (1)  it  is  necessary  to  recognize  the 
stochastic  effects  and  account  for  them,  C2)  instrument 
errors  are  important,  and  (3)  automation  of  the  model  struc¬ 
ture  development  stage  allows  quick  comparison  of  thousands  of 
both  simple  and  complex  models.  The  in-plane  moment  and  force 
coefficients  for  the  simulation  data  (without  boundary  layer 
transition)  are  explained  adequately  by  just  six  terms  and 
two  terms,  respectively.  The  model  for  the  ballistic  vehicle 
flight  data  during  boundary  layer  transition  is  more  complex, 
but  the  anomalous  nature  of  the  out-of-plane  moment  is  ex¬ 
plained  by  a  wind- fixed  flow  remaining  body  fixed  for  about 
0.2  sec  during  BLT.  Better  models  may  result  from  a  more 
detailed  simulation. 

6.2  RECOMMENDATIONS 

Though  some  significant  conclusions  may  be  drawn  from  this 
limited  application  of  the  techniques  to  one  of  the  BRV  flights, 
development  of  reliable,  better  detailed  models  will  necessarily 
require  extensive  application  of  these  techniques  to  a  wide 
spectrum  of  ballistic  vehicle  flight  data.  Data  from  a  single 
flight  just  does  not  contain  the  extent  of  information  required 
for  a  detailed  modeling,  partly  because  good  a  priori  information 
is  unavailable.  If  the  data  is  collected  on  several  vehicles, 
it  should  be  possible  to  more  accurately  correlate  BLT  mechanism 
with  vehicle  physical  characteristics.  This  would  provide  an 
effective  guide  for  future  designs  of  ballistic  reentry  vehicles 
which  would  minimize  impact  dispersion  resulting  from  the  BLT 
phenomenon. 

The  significant  modeling  aspects  of  ballistic  vehicles  can 
be  further  improved  by  more  detailed  analysis  of  instrumentation 
requirements.  Specialized  techniques  for  such  instrumentation 
specification  have  been  developed  [31]  and  can  usefully  be  applied 
to  such  ballistic  vehicle  requirements. 
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APPENDIX  A 

fstimation  of  windward  meridian  angle 

FROM  ONBOARD  SENSORS  USING  POLAR  COORDINATES 


Estimation  of  the  location  of  the  vehicle  velocity  vector  ts 
one  of  the  first  necessary  tasks  in  reconstructing  vehicle  states 
from  flight  data.  In  the  angular  rate  method,  it  is  assumed 
that  the  center  of  the  roughly  circular  path  described  by  the 
reentry  vehicle  (RV)  nose  coincides  with  the  velocity  vec  or. 

TO  determine  the  path  of  the  nose,  the  rate  gyro  measurements 

are  integrated,  using  zero  Initial  conditions,  to 

vehicle  pitch  and  yaw  Euler  angles.  The  mean  values  of  the 

the  nominal  velocity  vector. 

A.n  alternate  procedure  for  determining  the  velocity 
orientation  and  windward  meridian  angle  is  proposed  here,  lelocity 
orientation  is  determined  from  measurements  at  a  single  time  poin  . 
This  eliminates  the  need  for  rate  measurements  sufficient  to  a  o 
integration  over  one  cycle. 

The  procedure  is  based  on  two  assumptions: 


::ijfi?y^irc5of^rrentrt‘!irsft”irfy'^:ari"^ngie 
linearity  assumptions. 


Vehicle  acceleration,  as  measured  onboar 

accelerometers,  lies  in  the  wind  plane. 


Assumption  (1)  contributes  very  little  error  or  'h'  ^ 

of  attack  normally  encountered  «3"). 

,„t  appears  to  be  substantiated  j 

in  regions  o£  boundary  layer  transition  (BLT) .  Theie  , 
procedure,  if  it  is  to  be  used,  should  be  applied  to  flight  a 

prior  to  BLT. 


1 

If 


i 

I 
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Referring  to  Figure  A.l,  which  shows  a  rear  view  of  the  RV 

looking  along  the  velocity  vector,  the  following  definitions  are 
made : 


Figure  A.l  Trajectory  of  Vehicle  Nose 


The  vehicle  translational  velocity  is  assumed  normal  to  the  pitch- 
yaw  plane  and  passes  through  the  origin.  Initial  roll  angle  is 
assumed  to  be  zero. 

Five  unit  vectors  are  defined: 

A 

i^  -  tangent  to  the  path  of  the  nose 

A 

i,  -  traverse  to  the  vector  from  the  origin  to  nose 
^c  position 

A 

ij.  •  radial  along  vector  from  origin  to  nose  position 

A 

ig  -  in  direction  of  translation  due  to  pitch  (9)  rotation 

A 

i,  -  in  direction  of  translation  due  to  yaw  [ip]  rotation 


Two  angles  are  also  defined; 

ip  -  classical  precession  angle  relating  wind  plane  OA 
^  to  quasi-inertial  axes  (ig,  i^) 

a  -  total  angle  of  attack  between  vehicle  longitudinal 
axis  and  velocity  vector 

Based  on  Figure  A.l,  the  following  relationships  are  evident: 


V  =  vi^  =  ^i  +  9i^ 


(A.l) 


a>i<  i  I  oti- 

^c  r 


(A. 2) 


where 


V  is  the  velocity  of  the  nose  (normalized  by.  vehicle  length) 
The  transformations  between  polar  coordinates  i,  ,  i_ 

A  A  ^  C 

and  rectangular  coordinates  i^,  ig  are 


i^.^  «  cosn^^i^  -  sinij^^ig 


(A..3) 


ir  =  cosU^^ig  +  sini|;,.ig 


(A. 4) 


Combining  Eqs.  (A.#)  and  (A. 4)  into  A. 2),  letting  c  =  cos 
and  s  =  sin. 

A  •  ^ 

V  =  +  [-aij;|.si>^  +  aC'J;j.]ig  (A.  5) 

Comparing  Eqs.  (A.l)  and  (A. 5). 


ip  »  +  asrp^ 


(A.  6) 


(A. 7) 


Using  small  angle  assumptions 


^  =  qs(j)  +  rC(j) 


(A. 8) 


0  =  qc(t)  -  rs(J) 


(A. 9) 


where  the  vehicle  angular  velocities  are 


q  =  pitch  rate,  r  =  yaw  rate. 


and  where  the  vehicle  roll  angle  is  defined  in  Figure  A.  2. 


“  Y.q 


Figure  A. 2'  Vehicle  Roll  and  Windward  Meridian  Orientation 


Again,  it  is' assumed  that  the  vehicle  acceleration,  a,  lies  in 
the  wind  plane  OA.  Therefore, 


(A. 10] 


where  a„  and  a  are  measurements. 

y  z 


From  Figure  2,  the  wind¬ 


ward  meridian  angle  is 


=  180®  +  tan’-^  ^y/^z 


CA.ll) 


and  the  precession  angle  is 


The  roll  angle  (J)  is  easily  determined  from  roll  rate  measure¬ 
ments  except  for  an  arbitrary  initial  value. 

As  envisioned  in  the  application  of  this  procedure,  the  pitch- 
yaw  axes  are  assumed  initially  coincident  with  the  vehicle  y-z 
axes . 

Since  <j>  is  determined  directly  from  onboard  measurements 
and  <()  is  assumed  zero  (or  some  arbitrary  value)  ,  the  precession 
angle  can  be  directly  calculated  using  Eq.  (A.  12). 

Using  the  assumed  value  of  (j),  the  Euler  rates  (J ,  and 
(Ji  can  be  calculated  from  Eqs.  (A. 8)  and  (A. 9).  It  remains  to 
determine  a  and  a. 

If  Eqs.  (A. 6)  and  (A. 7)  are  squared  and  added,  Eq.  (A. 13)  results. 


+  9' 


2 


+ 


•  2 

a 


(A. 13) 


But , 


+  9^  =  (qC(J)  -  rsfj))^  +  (qs(l)  +  rc0) 

.  ,2  .  (A-l'') 

where  q  and  r  are  measured  quantities.  Therefore, 

7  2  2  •  2  *2 

q^  +  r^  =  +  a  (A. 15) 

Now,  ij.  in  Eq.  (-A.  15)  is 


v'c  V 


P  - 


(A. 16) 
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where  *  can  be  determined  numerically  from  successive  measure- 
w 

ments  of  vehicle  accelerations  using  Eq.  (A. 11),  and  p  is  the 
measured  roll  rate. 

Another  equation  in  addition  to  Eq.  (A. 15)  is  necessary  to 
find  a  and  a.  Multiplying  Eq.  (A.  6)  by  si|;^  and  Eq.  (A.  7)  by 
cip^,  and  adding  yields 

+  Qc'iJ  *  (A. 17) 

c  c  „ 


which  can  be  used  in  Eq. •  (A. 15)  to  yield 


01  *  /  ^ 

2  •  2 
+  r  -  a 

(A. 18) 

The 

pitch  and  yaw 

angles  can  now  be  determined 

9  =  aCii)^ 

(A. 19) 

\p  *  asip^ 

(A. 20) 

The 

•initial  vehicle  velocities  can  also  be  estimated  by 

u  =  V 

(A. 21) 

V  *  -Vij; 

CA.22) 

w  *  V9 

(A. 23) 
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APPENDIX  B 


MEASUREMENT  EQUATIONS 

This  appendix  gives  the  equations  for  measurements,  taken 
during  a  ballistic  reentry  vehicle  flight,  in  terms  of  the  equa¬ 
tions  of  motion  state  variables,  instrument  errors,  and  other 
related  terms. 

Equations  for  six  on-board  instrument  outputs  are: 

(1)  Axial  Accelerometer 

a^  -  u  -  vr  +  wq+(q2  +  ^^^^cg 

+  (pq  -  r)y^,g  +  (pr  +  q)Zcg  ^a^ 

(2)  Lateral  Accelerometer 

a^  *  V  -  wp  +  ur  +  Cpq  +  r)x^g  -  (r^^  +  P  -^cg 

+  (rq  -  p)z^„  +  b^ 

^  r.-  eg  a^^ 

(3)  Vertical  Accelerometer 


ay  -  w  -  uq  +  vp  +  (pr  -  q)Xj.g  +  (rp  +  p)y^^, 

-  ♦  p^)^cg  *  \ 

(4)  Roll  Rate 

p  ■  p  +  b 
‘^m  ^  p 

(5)  Pitch  Rate 
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